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Abstract 

In this thesis, we are interested by the stochastic control problems, we 
establish existence of an optimal stochastic relaxecl control for stochastic 
differential équation driven by a G —Brownian motion. 

We prove the regularity conditions realized by this control. 

Key words : 

G—brownian motion, coast function, relaxed optimal control. 
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Résumé 

Dans cette thèse on s’intéresse aux problèmes de contrôle stochas¬ 
tique, on établit l’existence d’un contrôle optimal relaxé pour une équa¬ 
tion différentielle stochastique gouvernée par un G—mouvement brow¬ 
nien. On démontre les conditions de régularité vérifiées par ce contrôle. 

Mots clés 

G—mouvement brownien, fonction coût, contrôle optimal relaxé. 


6 


Remerciements 


Je tiens à remercier tout particulièrement mon directeur de thèse le 
Professeur Hacène Boutabia pour l’honneur qu’il me fait en m’enca¬ 
drant et la confiance qu ’il m ’a accordée en acceptant de diriger ce travail. 
L ’orientation du sujet vers la théorie de contrôle était son idée, je le re¬ 
mercie pour ses conseils et pour tout le temps qu’il a consacré à cette 
thèse. J’aimerai également lui dire à quel point j’ai apprécié sa grande 
disponibilité et sa rigueur J’ai été extrêmement sensible à ses qualités 
humaines d’écoute et de compréhension. 

J’exprime mes plus sincères remerciements au Professeur Azzedine 
Benchettah, Professeur à l’université Badji Mokhtar de Annaba, et je 
lui assure ma profonde reconaissance pour avoir accepté de présider mon 
jury. 

Je souhaite exprimer ma gratitude et ma haute considération au Pro¬ 
fesseur Brahim Mezerdi, Professeur à l’université Mohamed Khider 
de Biskra, de m’avoir honoré par sa présence dans ce jury .Je le remer¬ 
cie de m’avoir rendue attachée à ce domaine de recherche depuis son 
encadrement de mon travail en magister. Mes premiers pas en calcul 
stochastique étaient sous sa direction et lors de toutes mes visites à Bis¬ 
kra il m’a écouté , il m’a donné de son temps.malgré ses responsabilités, 
mon invitation à V école d’hiver de probabilités qu’il a organisée à Biskra 
en janvier 2017 m’a fait un grand plaisir.et m’a aidé à améliorer mes 
connaissances . 

Je suis très reconnaissante aux Professeurs KANDOUCI Abdeld- 
jabbar et GUENDOUZI Toufik , Professeurs à l’université de Moulay 
Tahar de S aida d’avoir accepté d’être examinateurs de ma thèse. Je les 
remercie pour les efforts qu’ils ont fait pour faire partie de mon jury 
malgré la distance. 

Je remercie vivement le Docteur HADJI Mohamed Lakhdar, de 
l’université de Annaba pour avoir accepté d’être examinateur de cette 
thèse et pour l’intérêt qu’il a porté à ce travail, j"apprécie son attache¬ 
ment remarquable à la recherche. 

Au Professeur Boualem Djehiche, de l’institut royal de technolo¬ 
gie de Stockholm, Suède, dont je ne pourrai jamais le remercier assez 
pour tout son support pour réaliser ce travail, son écoute, sa disponibilité 
malgré ses tâches lourdes et ses explications valeureuses ; m’ont ouvert 
la porte pour atteindre mon objectif. 


7 


Grand merci à Salah Eddine Choutri, ça m’a fait plaisir qu’on 
partage ces résultats. 

Ma chère amie Dr. Goual Hafida, je te remercie pour ton soutien, 
en particulier ta précieuse aide à la rédaction de ce manuscrit. 

Ma grande amie Dr. Talhi Hamida, je trouve pas les mots pour 
te remercier. Ta présence à mes cotés et ton support durant toutes les 
phases de notre amitié sont inestimables... 

Ma chère amie Dr. Khodja Nawel, merci de tout coeur pour ton 
soutien illimité.. 

Mes vifs remerciements vont à toute ma famille pour tout ce qu’elle 
m’a donnée. 

Je remercie vivement tous mes amis (es) pour l’aide, le soutien et les 
encouragements. 

Je remercie également tous ceux qui, de près ou de loin, ont contribué 
à la réussite de ce travail. 


Dédicaces 

... Grâce à tes inestimables conseils, ta présence en chaque moment 
décisif de ma vie et tes encouragements je te rend hommage par ce mo¬ 
deste travail, à ta mémoire mon cher et regretté Papa je dédie cette 
réussite. 

... Toi, la lumière de ma vie, tu es toujours là pour nous, tes nom¬ 
breux sacrifices ont fait de moi ce que je suis, pour toi ma chère maman. 

... à mon mari Kamel, ta présence à coté de moi m’a éclairé la vie. 

... à mon unique soeur Meriam. 

... à mes frères Tarek et Mohamed. 

... à la mémoire de mon cousin regretté Ammar. 


9 


Introduction 

L’intérêt des problèmes de contrôle réside dans l’optimisation d’un 
certain critère de performance appelé fonction coût à l’aide du contrôle 
optimal et dans l’établissement des conditions nécessaires satisfaites par 
ce contrôle après avoir assuré son existence. En effet ; il existe de nom¬ 
breux domaines typiques tels que les conditions météorologiques, le cli¬ 
mat et les marchés financiers, où l’information est soumise à l’incertitude. 
Par exemple dans les problèmes de choix de portefeuille optimal en fi¬ 
nance, où les processus de volatilité et de prime de risque sont inconnus 
et difficiles à estimer à partir de données fiables, nous devons considérer 
une famille de différents modèles ou scénarios. Historiquement l’étude 
initiale du problème de contrôle est due à Lagrange qui a fondé le cal¬ 
cul des variations. Par la suite plusieurs auteurs ont mis l’accent sur 
les problèmes de contrôle optimaux parmi eux : Bismut, Benssoussan, 
Haussmann et Kushner. En 2006, Bahlali, Mezerdi et Djehiche [1] ont 
étudié le problème de contrôle optimal relaxé et ont démontré, sous cer¬ 
taines conditions de régularité, que chaque processus de diffusion associé 
au contrôle relaxé est une limite forte d’une suite de processus de diffu¬ 
sion associés aux contrôles stricts. 

Dans les problèmes de contrôle déterministe, la dynamique du sys¬ 
tème est modélisée par une équation différentielle ordinaire (EDO), par 
contre dans le cas stochastique l’évolution est décrite par un processus 
de diffusion solution d’une équation différentielle stochastique (EDS). 
D’un point de vue mathématique, on constate qu’il y a deux approches 
dans la résolution du problème de contrôle : la première approche est 
connue sous le nom de la programmation dynamique ou encore principe 
de Hamilton-Jacobi-Bellmann et la deuxième est le principe du maxi¬ 
mum appelé principe de Pontriagin. 

Parallèlement à cela, la théorie d’espérance non linéaire a été dévelop¬ 
pée et elle a eu beaucoup d’attention dans certains domaines d’applica¬ 
tion tels que la finance, la mesure des risques et le contrôle. Un exemple 
typique d’espérance non linéaire, appelé G—espérance a été introduit par 
Peng [34, 35]. La notion de la G—espérance s’est développée très récem¬ 
ment par Peng [34, 35] et a ouvert la voix à l’introduction de variables 
aléatoires G—normales, du G—mouvement Brownien et plus générale¬ 
ment des G—intégrales stochastiques de type Itô. Le G—mouvement 
Brownien est un processus stochastique avec des incréments station¬ 
naires et indépendants par rapport au passé et son processus de variation 
quadratique est, contrairement au cas classique, un processus non déter¬ 
ministe. En 2011, Soner,Touzi et Zhang [40], ont développé la G—analyse 
stochastique en introduisant la notion d’agrégation et la filtration uni¬ 
verselle. 
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Par la suite, Biagini, Oksendal et Paczka [6] ont étudié en 2014 le principe 
du maximum partant d’une équation différentielle stochastique contrôlée 
gouvernée par un G—mouvement brownien, sous l’hypothèse de l’exis¬ 
tence d’un contrôle stochastique optimal strict 

L’objectif de cette thèse est d’étudier le problème de contrôle 
stochastique partant d’une équation différentielle gouvernée par un G— 
mouvement brownien, où l’on montre en particulier l’existence d’un 
contrôle optimal relaxé vérifiant toutes les conditions de régularité né¬ 
cessaires pour garantir l’optimalité espérée. On s’intéresse à contrôler 
les systèmes soumis à l’incertitude ou à l’ambiguïté du modèle en raison 
d’informations incomplètes ou inexactes, ou de concepts vagues. Cette 
thèse s’articule autour de quatre chapitres successifs, dont voici la des¬ 
cription : 

Dans le chapitrel, on rappelle d’une part le cadre de la G—espérance 
en donnant les notions de base sur le G—mouvement brownien et d’autre 
part les notions de distributions et d’indépendance de variables aléatoires 
dans le cadre du G—calcul stochastique. 

Dans le deuxième chapitre, on fait un rappel du problème général du 
contrôle stochastique classique, en étudiant l’existence d’un contrôle op¬ 
timal strict minimisant un certain critère de performance appelé fonction 
coût. Dans certains cas le problème de contrôle ne possède pas une solu¬ 
tion optimale, en injectant par la mesure de Dirac, l’espace des contrôles 
stricts dans un espace plus large qui possède de bonnes propriétés to¬ 
pologiques notamment de compacité, on assure l’existence d’un contrôle 
optimal défini par une mesure de probabilité. 

Dans le troisième chapitre, on présente un aperçu général sur le tra¬ 
vail de Biagini et al [6], où ils ont vérifié le principe du maximum partant 
de l’hypothèse de l’existence d’un contrôle optimal strict associé à une 
G—équation différentielle stochastique. 

Enfin, on donne dans le quatrième chapitre nos principaux résultats 
sur l’existence d’un contrôle optimal relaxé. L’outil essentiel de notre 
contribution est le G—Chattering lemme, qui nous a permis de générali¬ 
ser le célèbre Chattering lemme classique au cas d’un espace d’espérance 
sous linéaire. On commence par définir l’espace des contrôles stochas¬ 
tiques relaxés à l’aide des outils récents du G—Calcul stochastique et 
on montre, avec détails sur les conditions de régularité, l’existence d’un 
contrôle optimal relaxé. Comme dans le cas classique, on montre que 
chaque contrôle relaxé est une limite d’une suite de contrôles stricts bien 
définie partant d’une G—équation différentielle stochastique. En fait on 
montre, à partir de l’approximation des trajectoires, que chaque proces¬ 
sus de diffusion lié au contrôle relaxé est une limite au sens fort d’une 
suite de processus de diffusion associés aux contrôles stricts. 
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Le présent travail de recherche a fait l’objet d’une publication inter¬ 
nationale intitulée : 

On Relaxed Stochastic Optimal Control for Stochastic Differential 
Equations Driven by G-Brownian Motion publiée à : ALEA - Latin Ame¬ 
rican Journal of Probability and Mathematical Statistics[39]. 


Par ailleurs, plusieurs communications en théorie de contrôle ont 
été faites : 

Communications internationales 

1— Approche variationnelle du principe stochastique de Pontriagin, 
Colloque International sur les Mathématiques Appliquées (CIMA’10, 7— 

9 novembre 2010, Guelma, Algérie). 

2— Approche variationnelle du principe stochastique de Pontriagin, 

Ecole CIMPA : Statistical methods and applications in Actuarial Science 
and Finance (CIMPA, 8 — 20 avril 2013, Marrakesh, Maroc). 

3— L’approche variationnelle du principe du maximum stochastique, 
Atelier international" Approximations dans les modèles stochastiques" 

(7 mai 2013, Béjaia, Algérie). 

4— Remarques sur le problème de contrôle stochastique, International 
Conférence of the Euro-Maghreb Laboratory of Mathematics and their 
Interactions (LEP21 — 2016, 27 avril-1 mai 2016, Hammamet, Tunisie). 

5— The G—Stochastic differential équations, Third International Confé¬ 
rence on Science, Health and Medicine 2016 (IRCSHM 2016,18 — 19 mai, 
Dubai 2016). 

6— The optimization in stochastic control theory, Third International 
Conférence on Science, Health and Medicine2016 (IRCSHM 2016, 18 — 19 
mai, Dubai 2016). 

7— The large déviations principle in G—Stochastic Calculus, Interna¬ 
tional conférence on stochastic Analysis, Stochastic Control 

and Application ,24 — 27 octobre, 2017, Hammamet, Tunisie). 

Communications nationales 

1— Le principe stochastique de Pontriagin, Congrès des mathémati¬ 
ciens algériens (CMA—2012, 7 — 8 mars, Annaba, Algérie). 

2— Conditions nécessaires d’optimalité dans le contrôle stochastique, 
2 îemes Journées sur les problèmes inverses : Théorie et applications (JIP’13,28— 
30 octobre, Annaba, Algérie) 
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3— L’optimalité dans la théorie de contrôle stochastique, Conférence 
sur les Problèmes Mathématiques non linéaires (ProMa’2014, 11 — 13 
mars 2014, Annaba, Algérie). 

4— Propriétés d’un contrôle stochastique optimal, Congrès des ma¬ 
thématiciens algériens(CMA-2014, Tlemcen, Algérie). 

5— L’approche variationnelle en théorie de contrôle stochastique, Jour¬ 
nées Jeunes Chercheurs (J JC, 30 septembre-1 octobre 

2014, Annaba, Algérie). 

6— The optimality in Relaxed Stochastic Control Problems, Premier 
Séminaire National de Mathématiques(SNM’01,13 décembre 2016, Constan- 
tine, Algérie). 


13 


Chapitre 01 

Généralités sur le G— calcul 
stochastique 


14 


1 Généralités sur le G—calcul stochastique 

L’objectif de ce chapitre est de rappeler quelques notions de base 
qui seront utilisées tout au long de ce travail. 

Dans ce chapitre, on introduit quelques notations et préliminaires, 
notamment sur la théorie d’espérance sous linéaire, le G—mouvement 
Brownien hé à la distribution G—normale. Pour plus de détails, le lecteur 
pourra se référer à [23,27,33, 36, 37]. 

1.1 Espérance sous linéaire 

Soit Q un ensemble donné et soit Ti un espace linéaire de fonc¬ 
tions à valeurs réelles définies sur D, tel que c G PL pour chaque constante 
c et |X| G PL si X G PL. PL est considéré comme l’espace des variables 
aléatoires. 

Définition 1.1 Une espérance sous linéaire sur H est une fonction 
E : H —> R satisfaisant les propriétés suivantes : Pour tout X , Y G H, 
on a 

1 Monotonicité : Ê [X] > Ê [Y] si X > Y. 

2 Préservation des constantes : E [c] = c, pour tout c G M. 

3 Sous-additivité : Ê[X\-Ê [Y] < Ê [X - Y}. 

4 Homogénéité positive : E [XX] = À E [X], pour tout À > 0. 

le triplet (çi,PL,È^j est appelé espace d’espérance sous linéaire. 

Si seulement (1) et (2) sont satisfaites, E est appelée espérance non 
linéaire. 

Si l’inégalité (3) est une égalité, alors E est une espérance linéaire 
classique, i.e., E est une fonction linéaire satisfaisant (1) et (2). 

Remarque : En fait (3) et (4) impliquent la propriété de convexité 
suivante : 

Ê [aX + (1 — a) F] < aÊ [X] + (1 — a) Ê [F] pour a G [0,1]. 

Notons que la propriété (4) est équivalente à la propriété suivante : 

Ê[\X) = \ + È[X ] + \~Ê[-X] pour A G M, 

où À + = max (À, 0) et À - = max (—À, 0). 

Dans toute la suite on désigne par C^Lip (M 71 ) l’espace des fonctions 
réelles définies sur K n telles que 

|(p (x) — (p (y)| < G ^1 + |x| fc + |y| fc ^ |x — y\ pour tout x, y G 
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où k G N* et G > 0 dépendant uniquement de <p. 

Théorème 1.2 ([7,34]) : Soit E une espérance sous linéaire défi¬ 
nie sur un espace linéaire H , alors il existe une famille de mesures de 
probabilités V sur (0,0(0)) telle que 

Ê [X] = sup Ep [X] pour X G T~L, 

PeP 

où B (O) désigne la tribu borélienne sur O. 

Naturellement, on peut définir une capacité de Choquet sur T-L par 

c (A) := sup P ( A) , A G B (O), 

PeP 

Pour plus de détails, le lecteur pourra consulter [10, 21,43] . 

Définition 1.3 On dira qu’un ensemble A E B (O) est polaire si et 
seulement si c(A) = 0 et qu’une propriété a lieu quasi-sûrement (q.s. 
en abrégé) si elle a lieu en dehors d’un ensemble polaire. 

Remarque : Une propriété est vraie q.s. si elle est vraie presque 
sûrement pour chaque P G V. 

Grâce à la représentation de la G—espérance et à la notion "quasi- 
sûre" introduite, la théorie des processus stochastiques en temps continu 
dans le cadre de la G—espérance s’est développée, en particulier la for¬ 
mule d’Itô, certaines inégalités stochastiques et les équations différen¬ 
tielles stochastiques gouvernées par un G—mouvement Brownien (G—EDS 
en abrégé) peuvent être établies dans le sens "quasi-sûre". 
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1.2 Distributions et indépendance 

Parallèlement aux concepts du cadre classique, Peng [35] a établi 
les notions de distributions et d’indépendance pour les variables aléa¬ 
toires dans ce nouveau contexte. Néanmoins, ces notions sont moins 
probabilistes mais plutôt fonctionnelles, et elles s’expriment à l’aide des 
familles de fonctions tests de l’espace Cij jVp (R”), n > 1. 

Définition 1.4 Soit X = (Xi,...,X n ) un vecteur aléatoire, où Xi G 
H, i = 1 La distribution de X est donnée par la fonctionnelle 

Fx [•] suivante : pour chaque p G C^Lip (R™), 

F X M :=Ê[<p(X)]. 

Le triplet (R n , C^Lip (R n ), F^) forme un espace d’espérance non li¬ 
néaire. 

De plus, si X' est un autre vecteur aléatoire n—dimensionnel et pour 
chaque p G C'/.Lip ( R n ), Fx[p\ = Fx> [p], X et X' sont dits identique¬ 
ment distribués. 

Remarques : 

- On peut prouver (voir [7,30],) qu’il existe une famille de mesures 
de probabilité {F[^(.)} ee0 définie sur (R n ,£>(R n )) telle que 



Ainsi Fx [p] caractérise l’incertitude de la distribution de X. 
Notons que si la distribution F^ de X G LL n’est pas une espérance 
linéaire, alors X a une distribution incertaine. La distribution de X a 
les quatre paramètres typiques suivants : 

JI = Ê[X], p = -Ê [-X ], if 2 = Ê [X 2 ] et g 2 = -Ê [-X 2 ] . 

Les intervalles [jf l 1 ] et [F 2 , a 2 ] caractérisent la moyenne incertaine 
et la variance incertaine de X respectivement. 

Proposition 1.5 : Soient X,Y G H telles que Ê[Y] = —Ê[—Y) ( 
i.e., Y n’a pas une moyenne incertaine). Alors, on a : 

Ê[X + aY]=Ê [X] + aÊ [Y]. 

En particulier, si Ê [F] = — Ê [—F] = 0, alors Ê [X + a F] = Ê [X], 
Preuve : On a 

Ê [aY] = a + Ê [F] + oTÊ [-F] = a + È [F] - oTÊ [F] = aÊ [F] 
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pour a£l. Ainsi 


Ê [X + aY]<Ê [X] + Ê [aY] = Ê [X] + otÊ [F] = Ê [X] - Ê [-«F] 

< Ê [X + aY] . 

■ 

Proposition 1.6 : Soient <p G Ci ! i 4 P (W n+n ) et Y G H n un vecteur 

aléatoire défini sur un espace d’espérance sous linéaire PL, Ê'j . Alors 
on a : 

1) Pour tout x G l’application <p (x ,.) G C^up (®P)- 

2) L’application Ê [p (., F)] G C^u p (K m ). 

Preuve : 1) En effet, pour chaque (x, y ), (u, z ) G M m x 1", on a 

I <p(x,y) ~ <p(u,z)\ <C \(x,y) - (■ u,z)\ (l + \(x,y)\ k + |(u, 2 )| fe ) , 

C>0,k G N, 

et ainsi 

I <p(x,y) - p(x,z)\<C\y - z\ (l + |(x,y)| fc + |(x, 

<C\y - z\ ^1 + (|x| 2 + |y| 2 ) 2 + (|x| 2 + |-| 2 ) 2 ^ 
<C\y-z\ (l + 2* (\x\ k + \y\ k ^j +2 2 (jx| fc + 

<C'\y-z\ (l+ (|y| fc + M fc )) , 

où on a utilisé l’inégalité (a + b) p < 2 P (a p + bP) pour a, b > 0 et C' — 

C max ^2^, 1 + 2^ +1 \x\ l ^j . Par conséquent, (p ( x ,.) G Cip p (M n ). 

2) De même, on a 

|(p (x, Y) - p (u,Y )| <C\x — u\ (l + 2* (\x\ k + \Y\ k ^j + 2^ (\u\ k + |F| fc ^ , 


ce qui implique que 


Ê[p (x,Y)]-Ê[p(u,Y)] 


< Ê \ p(x, F) -p(u,Y )| 


<C\x-u\ (l + 2* (\x\ k + Ê (|F| fc )) +2^ {\u\ k + Ê (|F| fc ))) 
< C” \x — u\ ^1 + \x\ k + , 
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où C" = C max ^22,1 + 22 +1 E (jF| fc ^, ce qui signifie que Ê [<p (., F)] 

appartient à Ci t u p (R m ). 

■ 

Définition 1.7 Dans un espace d’espérance sous linéaire (^d,TL, Ê^j, 
un vecteur aléatoire Y G H n est dit indépendant d’un autre vecteur 
aléatoire X G H m sous Ë si 


Ê[<p(X,Y)] 


E 


Ê [p (x, F)] |x=A' 


y p g q Mp (R m+n ). 


Définition 1.8 Soient X\ et X 2 deux vecteurs aléatoires à n—dimension 
définis respectivement sur deux espaces d’espérance sous linéaire TL\, Ê\ 

et ^ 2 ; H 2 , Ê 2 ^j • On dira que X 1 et X 2 sont identiquement distribués et 

on notera X 1 = X 2 , si 

Ê 1 [p (Xi)] = È 2 [p (X 2 )] pour tout p G C iMv (R"). 


Si X est indépendant de X et X 
copie de X. 


d 


X, alors on dira que X est une 


Remarque : Dans un espace d’espérance sous linéaire la condition 
"Y est indépendant de X " ne signifie pas automatiquement que "X est 
indépendant de Y ([43]). 
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1.3 Distribution G—normale et G—mouvement Brow¬ 

nien 

Après la définition de base ci-dessus on introduit maintenant la no¬ 
tion de la distribution G—normale. 

Définition 1.9 Un vecteur aléatoire à d—dimension X = (Xi, X d ) 
dans un espace d’espérance sous linéaire LL, Ê^j est dit G—normalement 
distribué si pour tout a, b > 0 : 

aX + bX = Va 2 + b 2 X, 
où X est une copie indépendante de X, et 

G (A) := l -Ê[(AX,X)} : § d ^ K, 

ici Sd désigne l’ensemble des matrices symétriques d’ordre d. 

Dans [7,15,37], Peng montre que X — (Ad, ...,X d ) est G—normalement 
distribué si et seulement si u (t, x ) := Ê \<p (x + VtX)] , (t, x ) G [0, oo) x 
W l , p G Ci^ (K d ) est l’unique solution de viscosité de l’équation para¬ 
bolique aux dérivées partielles, appelée G—équation de la chaleur, sui¬ 
vante 

f d t u ( t , x) = G ( Du ( t , x )), ( t , x) G [0, oo) xK d , . 

( u (0, x) = p (x) ’ 

où Du (t, x) est la matrice Hessienne de u (t, x). 

G (.) : S d —» R est une fonction monotone et sous linéaire sur Sd, à 
partir de laquelle on peut déduire qu’il existe un sous ensemble E G §} 
fermé, borné et convexe tel que 

G (A) = - sup tr (AB). 

2 BeT. 

On écrit X ~ J\f (0; S). 

Remarque : Le cas réel (d = 1) correspond à S = [<t 2 ;< J 2 ] et G = 
Ga,cj étant la fonction sous linéaire paramétrée par a et a : 

G (a) = - (<r 2 a + — a 2 a ~) , a G I, 

oùW 2 = Ê [X 2 ] et a 2 = —Ê [—X 2 ]. Dans ce cas on écrit X ~ J\f (0; [çr 2 ; cr 2 ]). 

Corollaire 1.10 Dans le cas où a 2 = a 2 > 0, AT(0; [a 2 ] o 2 ]) est juste 
la distribution normale classique J\f(0',a 2 ). 
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Preuve : En fait, la solution de l’équation parabolique aux dérivées 
partielles d t u = G (d xx u) , (t, x ) G [0, oo) x K avec la condition de Cauchy 
u (0, x) = p ( x ) devient l’équation de la chaleur classique 


d t u = y dl x u, u | t=0 = (p, 


où la solution est 


u ( t , x) = 


1 


+oo 


VËïrcrH . 

Ainsi, pour chaque p, 

Ê[<p(X)] = u( 1,0) = 


ï ’ (!,)exp ( Jl wf) iy - 


+oo 


\J‘1'ïïo 2 


P (y) exp ( -y 2 ) dy 


Dans les deux situations suivantes le calcul de Ë [p (X)] est facile : 
•Pour tout p convexe, on a 


u ( t , x) 


+oo 




P (y) exp 



dy. 


En effet pour tout t > 0 fixé, remarquons que la fonction u ( t , x) 
Ë [(p (x + a/ÂA)] est convexe, puisque 


u (t, ax + (1 — a) y) = Ê p [ax + (1 — a) y + VtX^j 


< aE 


( 


p\x + 


+ (1 — a) Ê 


p[y + 


Vtx\ 


= au ( t , x) + (1 — a) u ( t , y) 


Il en résulte que (d xx u) = 0 et par conséquent la G—équation de la 
chaleur (1) devient 


d t u 



u h =0 = (p. 


•Pour tout (p concave, on a 


u ( t , x) 


+oo 


\[Ëkq} 


p (y)exp 



dy, 
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en particulier 

Ê [X] = Ê [-X] = 0, Ê [X 2 ] = â 2 , -Ê [-X 2 ] = a 2 , 
et 

Ê [X 4 ] - 6â 4 , - Ê [-X 4 ] = 6a 4 . 

On donne maintenant la définition du G—mouvement Brownien d— dimensionnel. 


Définition 1.11 Un processus ( B t ) t>0 à d—dimension défini sur 
est dit G—mouvement Brownien de dimension d (réel dans 
le cas où d — l) si les propriétés suivantes sont satisfaites : 

(■ i ) B 0 = 0 

(ri) pour chaque t,s > 0, l’accroissement B t+S —B t est N (0, s£) — distribué 
et est indépendant de (B tl ,..., B tn ), pour tout n G N et 0 < ti < 

... < t n < t. 

On notera 

Bf (a, B t ) pour tout a — (ai,..., ad) T G W 1 . 

Selon la définition ci-dessus, on a la proposition suivante qui est impor¬ 
tante dans les développements ultérieurs. 


Proposition 1.12 ( voir [17,38]) Soit (B t ) t>0 un G—mouvement 
Brownien à d—dimension défini sur un espace d’espérance sous linéaire 
(fi,, EL, Ê ^ . Alors (Bf) t>0 est un G a —mouvement Brownien réel de fonc¬ 
tion génératrice G a (a) = \ ((7^ t q + — gfi aaT a~) où (f 2 aaT = Ê [(a,Bi) 2 ] 
et a 2 aaT = —Ê [— (a, i?i) 2 ] . En particulier, pour tout t,s > 0, Bf +S — 
B t ~-A/"(0; [sa 2 aT ;sâ 2 aT ]) . 

Notons que d’après le cas unidimensionnel, on a pour toute fonction 
convexe ip : 


+oo 


ê [f {b; + , - b^] 


27rsa 2 


(p(y)exp ( - 0s _ 2 ) dy, 


et pour toute fonction concave <p et a 2 aaT > 0, on a 

+oo 


Ê [p (Bf +S - Bf)\ = 


1 


27rsa 2 


<P (d) ex P - 


2sa 2 


dy. 


-aa* — oo 


En particulier, d’après Peng [34,36,38], on a 

ê [(b? - b:ÿ] = âi aT (t — s) ,ê [{Bf - b:ÿ] = 

Ê [- ( 5 “ - B a s Ÿ] = -air (t-s),Ê[- ( 5 “ - Bf) 4 } 


3c7 Lr (t-S) 2 , 

= — 3a 4 a T (t - s f . 
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Chapitre 02 

Problème de.contrôle 
stochastique 
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2 Problème de contrôle stochastique 

Le problème de contrôle constitue une partie essentielle dans le 
domaine des mathématiques appliquées, il consiste à optimiser un critère 
de performance d’un système dynamique. Dans le cas déterministe, la 
théorie de contrôle modélise les problèmes d’évolution où la dynamique 
est décrite par une équation différentielle ordinaire, alors que le cas du 
contrôle stochastique traite la contrôlabilité des systèmes dynamiques 
évoluant dans des conditions d’incertitude. Dans ce cas l’évolution est 
modélisée par une équation différentielle stochastique (EDS d’Itô). 

2.1 Problème de contrôle déterministe 

On considère un système dynamique où l’évolution est guidée 
par l’équation différentielle ordinaire contrôlée suivante : 

{ Ht u 
Xq = X 0 

où / : M + x R p x R d —* M p et u : R + —» Aci 11 sont deux fonctions 
données. On considère aussi U C F (R + , A), l’ensemble des contrôles ad¬ 
missibles, où F(R + , A) désigne l’ensemble des fonctions de M + vers A et 
on suppose que pour tout u G U, il existe une solution x u du système 
ci dessus. Soit J une fonctionnelle de u définie sur U par J(u) — 
où $ : U — > R, T étant un temps terminal fixé dans R + , et on cherche 
un contrôle u* EU réalisant u* = inf ue u J(u), minimisant le coût J. 


Exemple : On considère un système dynamique, où l’évolution est 
modélisée par une équation linéaire : 


dx u 

dt 


a t x t + b t u t . 


Le côut est défini par : 

r T 

J(u ) = / m s x 2 s + n s u 2 s ds + dx 

J o 

avec mn s > 0, n s > 0, d > 0. On pose = (x“, y“) avec : 

Vt — f (m s Xg + n s u 2 )ds et $(X^) = + d(x^) 2 . 

J o 

Dans ce cas on a J(u) = $(Xÿ) (voir Mazliak[30]). 


24 


2.2 Problème de contrôle stochastique 

On se place dans un espace probabilisé filtré (O, F, (F t ) teR +, P), 
satisfaisant les conditions habituelles et on se donne B = (B t ) teR + 
un mouvement brownien d— dimensionnel et A un borélien de R. Dans 
les problèmes de contrôle stochastique, on distingue deux formulations 
comme suit : 

i) Si l’espace de probabilité filtré et le mouvement brownien sont donnés 
à l’avance et les contrôles sont adaptés à cette filtration, on a une for¬ 
mulation forte. 

ii) Si l’espace de probabilité filtré et le mouvement brownien sont des 
parties de contrôle on va obtenir une formulation faible. 

Les contrôles admissibles sont l’outil de base dans notre problème de 
contrôle. 

Définition 2.1 On appelle contrôle admissible tout processus 
leurs dans A, mesurable et (F t )tës+ —adapté. 

On rappelle qu’un processus u : D —* F([0, T], A) est ( F t ) teR + 
si Ut est une application F t — mesurable. 

Remarque : Le processus u doit être (F t ) — adapté pour assurer la 
connaissance de l’histoire du processus jusqu’à l’instant t. 

2.2.1 Position du problème 

Cas où le terme de diffusion est contrôlé 

Etant donné un système dynamique où l’état est modélisé par l’EDS 
gouvernée par un contrôle : 

dX t = b(t , X t , u t ) + cr(t, X t , u t )dB t 
Xo = x, t G ]0, T] , x G R m 

où b : [0, T] x R m x A —» M m et a : [0, T] x R m x A —* Af TOX d(R) sont 
deux fonctions boréliennes, dérivables et à dérivées continues et bornées, 
Af mX d(R) étant l’espace des matrices d’ordre (mxd) a coefficients réels. 
On suppose en outre que : 

-La fonctionnât, x, a) est continue en a et uniformément continue en 
t et en x. 

-Le processus X est contrôlé a chaque instant t G [0, T] a l’aide de la 
valeur du paramètre a prise dans le borélien A. 

-Les fonctions b et a sont à croissance linéaire : 

|6(t,x,a)| 2 + |<j(t, x, a)\ 2 < C(1 + |x| 2 ) (3) 


à va- 
—adapté 
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où C est une constante positive indépendante de a et de t. 

Notons que b et a sont dans ce cas lipschitziennes, il existe alors une 
constante positive K indépendante de t et a telle que : 

Vx, y G R m : | b(t, x, a) — b(t, y, a) | 2 + | a(t, x, a) — a(t, y , a) | 2 < K\x — y\ 2 

( 4 ) 

Sous les les hypothèses : (3), (4), l’EDS (2) admet une unique solution 
X, d’après le théorème d’existence et d’unicité. 

Notons que le processus X est progressivement mesurable autrement 
dit non-anticipatif. 


Définition 2.2 Le problème de contrôle est fondé sur un critère de 
performance appelé fonction coût définie par la forme suivante : 


rT 


J(u) = E[ l(t, X t , u t )dt + g(X T )\ 


( 5 ) 


où g : M m —^ K et l : [0, T] x R m x A —» R sont des fonctions réelles. 
Il est clair que la variable aléatoire Xt est associée à un contrôle admis¬ 
sible u pour lequel on définit le coût J{u). 

L’objectif visé par le problème de contrôle est de déterminer l’infimum 
pour J et par suite, trouver un contrôle s’il existe, qui le réalise pour 
obtenir l’optimalité espérée. On doit donc prouver d’abord l’existence 
d’un contrôle admissible u* vérifiant pour tout x G R m : 

V(x):= inf J u (x) = J u \x) 

ueu 

La fonction V est appelée coût optimal ou bien fonction valeur et le 
contrôle u* est dit optimal. 

Cas où le terme de diffusion n’est pas contrôlé 

Dans le cas où le terme de diffusion n’est pas contrôlé, l’état de notre 
système dynamique est décrit par l’EDS suivante : 


dX t = b(t , X t , u t )dt + a(t, X t )dB t , t G ]0, T] 
X 0 = x, x G R m 


( 6 ) 


où b, a sont deux fonctions boréliennes satisfaisant les conditions décrites 
auparavant. Notons que dans cette EDS on ne contrôle que le drift. La 
fonction coût est définie, dans ce cas, de la même manière que dans (5). 
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Exemple : Dans le problème d’arrêt optimal, l’objectif est de mini¬ 
miser la fonction coût : 

J(u,t)=E[ f l(t,X u u t )dt + g(X T )\, 

J o 

sur les contrôles admissibles u et les temps d’arrêt r. 

On considère l’exemple de vente d’un actif : une personne désire vendre 
un actif ou une ressource (maison,action,...), l’évolution du prix de cet 
actif est modélisée par l’EDS : 


dX t = r.X t dt + a.X t dB t . 

On suppose qu’à la vente de l’actif, il y a un coût de transaction fixe 
a > 0. Si la personne décide de vendre l’actif à la date t, le bénéfice net 
de vente sera exp(— ip t )(X t — a) où <p > 0 est le facteur d’inflation. Le 
problème est de trouver un temps d’arrêt qui maximise la valeur espérée 
du bénéfice net à savoir supE[exp — (p T \(X T — a) et de calculer le profit 

T 

espéré. 

Comparaison entre le contrôle déterministe et le contrôle 
stochastique 

Dans le contrôle stochastique, la condition d’adaptation du contrôle 
est nécessaire pour que l’intégrale stochastique soit bien définie. D’autre 
part on ne peut pas optimiser le coût J pour tout weîî comme dans le 
cas déterministe, c’est pour cette raison qu’on minimise un coût moyen. 

2.3 Existence d’un contrôle optimal ordinaire 

On considère le problème de contrôle déterministe suivant : A — 
{ — 1,1} représente l’ensemble des contrôles et U représente l’ensemble 
des fonctions continues par morceaux, nommées fonctions contrôles, u : 
[0,1] —>■ A. On se donne l’équation différentielle ordinaire suivante : 

f dx u [t) = u(t)dt, t E [0,1] 

\x(0) = Xq 

et on définit la fonction coût par J(u) = f 0 (x u (t)) 2 dt. 

Remarquons d’abord que inf u J(u) = 0. On considère une suite (u n ) ne pj 
telle que Vn G N* on a : 

U n {t) = (-l) fc , £ < t < ti, 0 < k < n - 1. 

n 
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Il est clair que \x Un (t)\ < —,Vt G [0, l],et que J(u ) < 

Si ou suppose que x u (t) = 0 pour tout t G [0,1] ou obtient u t = 0 pour 
tout t G [0,1], ce qui est impossible. On conclut qu’il n’existe pas un 
contrôle u G U réalisant J(u ) = 0. 

La suite ( u n ) n’admet donc pas de limite dans l’espace des contrôles U, 
tenant en compte que cette limite doit être le candidat naturel d’opti¬ 
malité, on en déduit que notre problème de contrôle n’a pas de solution 
optimale, c’est pour cette raison qu’on va définir un espace plus large 
vérifiant certaines conditions de régularité dans lequel cette limite existe. 

2.4 Problème de contrôle strict 

On considère l’EDS contrôlée (2) sur l’espace (O, F, (F t ) teR +, P), avec 
les mêmes hypothèses précédentes. Le générateur infinitésimal L associé 
à l’EDS (2) appliqué sur les fonctions :/ G Cf(M m ,M) est défini par : 


Lf(t,x,u ) = 


d 2 f f)f 


où cr T désigne la transposée de a. 


Définition 2.3 On appelle contrôle strict la donnée du septuplet 
a = (O, F, F t , P,u t ,Xt,x), où (u t ) est un processus à valeurs dans A, 
progressivement mesurable par rapport à (F t ) teR + et ( X t ) représente 
la solution de l’EDS (2) avec la valeur initiale x G R m , telle que le 
processus : 


f{xt) ~ f{x) 



Lf(s,X s ,u s )ds , 


est une P—martingale pour toute f G C 2 et pour tout t > 0. 


Remarque : Si le contrôle u t est constant, les conditions données 
sur b et a assurent que notre problème de martingale admet au moins 
une solution, ce qui implique une existence faible de la solution de l’EDS 
(2). Dans ce cas les contrôles associés sont dits contrôles faibles car l’es¬ 
pace de probabilité et le mouvement brownien peuvent être modifiés selon 
Ut- 

L’unicité trajectorielle a été étudiée par El Karoui et al [11], où ils 
montrent que les problèmes de contrôle fort et faible ont les mêmes 
fonctions valeurs. 
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2.5 Optimalité des contrôles adaptés à l’état du sys¬ 
tème 

Le but de cette étude est la minimisation du coût concernant les 
contrôles adaptés à l’état du système. On associe pour tout u G U, la 
filtration naturelle (F't 1 ') du processus X (Le. Ffi = a(X(s),s G [0, t[)) 

et on cherche à montrer que Ffi — F t . Par la suite, on utilise ce résultat 
pour démontrer que l’espace U* est dense dans U, où U* est donné par : 

U* = {u G U : u(t ) est Ffi— mesurable p.p.}. 

Il résulte alors que inf J(u) = inf J{u). 

ueu* uÇiU 

On a besoin de la définition suivante d’un bruit additif : 

Définition 2.4 On appelle bruit additif tout processus stochastique 
continu adapté noté br(t), à valeur dans K m ,tel que : 

sup \br(t)\ 2 < C. 

0<t<T 

Lemme 2.5 Soient u G U*et (F t ) la filtration naturelle de ( B t ). On 
considère l’équation : 

fx(t) = x 0 + /g b(s,X s ,u 3 )ds + br(t) 

\x 0 =.x 

alors on a F t = Ffi pour tout t. 

Preuve : Il est clair que Ffi C F t pour tout u G U. Il reste à dé¬ 
montrer que F t C Ffi. Pour cela, on définit 

p{t) := / (T(s, X s )dB s = X t - x - / b(s, X s , u s )ds, 

J o J o 

Le processus p(t) est une (F t ) —martingale continue et le processus 
f* aa T (s, X s )ds est bien défini. Il résulte que l’intégrale stochastique 
I{p) = So^~ 1 (s,X a )dfj,{s) a un sens. 

Le fait que les intégrales stochastiques sont des limites au sens de 
Riemann, cela nous ramène au résultat suivant :/(t) = B t . 

Maintenant on considère u G U*. D’après la définition de p(t), il est 
clair que (p(t)) est -adapté. En utilisant l’inclusion Ffi C F t , d’où 

l’adaptation de I{t ) et B t par rapport à X * ; Il résulte que F t C Ffi, 
d’où le résultat demandé. ■ 
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lemme 2.6 Sous les mêmes hypothèses du lemme précédent, l’espace 
U* est dense dans U. 


Preuve : Soient u G U et K = -F . On considère le processus à 
valeurs dans A défini par : 

^ çnK 

u k (t) = - u(s)ds si t G [nK, (n + 1 )K[ 

k J(n-l)K 

pour n — 1 , N — 1 . 

D’après les propriétés de A, on déduit que u k G U . De plus la suite 
(uk) converge vers u dans L 2 ([0,T]). 

Il restera déduire que u k EU*. Soient X k la trajectoire correspondante 
à Uk, F kt la tribu engendrée par X k (t) et 


Sk(t)= / a (s, Xk(s))dB s 


= X k (t) -x 0 - b(s , X k {s), u s )ds , 


et : 

B{t) = f* cr _1 (s,X fc(s) )d// fc (s). 

D’après la définition de Uk, on déduit que pLk(t) es t Fk,t ~mesurable 
pour t G [0, A;].Puisque F, = F k)t pour tout 0 < t < k, alors Uk(t) est 
Fk,t ~mesurable pour t G [k,2k\. D’après la formule ci-dessus et étape 
par étape, on aura F t = pour tout t , d’où le résultat désiré. ■ 

On alors le théorème suivant : 


Théorème 2.7 Sous le hypothèses précédentes on a 

inf J(u) = inf J(u). 

uEU uEU* 

Preuve : Il suffit de démontrer que J est continue au sens de L|.(0, T), 
ce qui revient à montrer que pour toute suite (uk) convergente vers u 
dans L 2 f (0, T) on a lim J(uk) = J(u) dans K. 

k—>oo 

On considère d’abord le cas d’un bruit additif : 

f X t = x 0 + /J b(s, X s , u s )ds + br(t ) 

[X 0 = x 

On a : 

X k {t) - X(t) = / b{s,X k {s),u k {s))ds - / b(s, X s , u s )ds 
J 0 J 0 
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et d’après la condition de Lipschitz de b , on aura 



où l’on a posé 



Par suite on obtient <Pk(T) —» 0. 


D’après la condition de croissance linéaire sur J, on a l’inégalité : 



\J(u k ) - J(u )| <C[E (\X k (t)\ + \u k {t )| + 1)(| X k (t) - X(t )| + | u k (t) - u(t)\)dt 


+E(\X k (T)\ + | X(T) + 1|| X k (T) - X(T)|)], 


d’où le résultat.■ 

Remarque : Si la fonction coût J atteint son minimum pour le 
contrôle u G U*, u est aussi un contrôle optimal dans U donc il doit 
vérifier les conditions nécessaires d’optimalité dansU. 

2.6 Quelques exemples sur le contrôle stochastique 

2.6.1 Contrôle de risque sensitifs 

Soit ip : K —ï K une fonction monotone, on définit le coût par : 




J{u) = Eip[ / l(t, X t , u t )dt + g(X T 


Notons que X t est la solution de l’EDS modélisant le système dynamique 
à étudier prise à l’instant t. Ce genre de problème est intéressant dans 
la théorie économique. D’habitude les fonctions de coût sont de la forme 
Lp e (x) = 9 exp (9x) 



ln x si : 9 = 0 


si : 9 ^ 0 


2.6.2 Contrôle ergodique 


On considère le problème de contrôle stochastique dans lequel l’ho¬ 
rizon est infini, le coût est défini par : 



Ce problème est nommé un problème ergodique. 
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2.6.3 Portefeuille de Merton 

On considère le problème du choix de portefeuille de Merton où l’ho¬ 
rizon est fini. On se place dans un marché financier à deux actifs, l’actif 
sans risque représente le compte d’épargne et l’autre représenté par une 
action est considéré risqué. Le processus de prix de l’actif sans risque 
évolue selon l’équation ds® = rs^.dt et celui de l’actif risqué évolue selon 
l’EDS ds t = fjLS t dt + a s t .dB t , où r, p, et a sont des constantes avec : 
a > 0 . 

Soit un agent qui investit dans ces deux actifs avec une proportion de sa 
richesse. A la date t, cette proportion est notée 7p dans l’actif risqué et 
(1 — 7r t ) dans l’actif sans risque. 

Son processus de richesse évolue selon : 

dX t = (1 - 7T t )^-ds° t + 7 Tt—dst 

s t s t 

= [ 7 T t X t p + (1 - 7r t )X t r]dt + 7 Y t X t adB t . 

Le contrôle est le processus 7 r sont à valeurs réelles. Le critère écono¬ 
mique consiste à maximiser l’espérance de l’utilité de la richesse termi¬ 
nale à un horizon fini T : sup^ E[u(Xt)\ où u est une fonction d’utilité, 
prenons par exemple u(x) = —, : p G]0,1]. 

2.7 Le problème de contrôle relaxé 

D’une façon générale, le problème d’optimisation en théorie de 
contrôle stochastique, ne possède pas de solution optimale en l’absence 
d’hypothèses supplémentaires de convexité sur la dérive b et le coeffi¬ 
cient de diffusion a. Pour cette raison on injecte l’espace des contrôles 
stricts U, dans un espace plus large qui possède de bonnes propriétés de 
compacité et de convexité et on pose le problème de contrôle optimal 
relaxé sur cet espace noté 3Î. 

Notons que K est muni de la topologie de la convergence stable (topologie 
faible). L’existence d’un contrôle optimal pour les problèmes modélisés 
par les EDS dont a ne dépend pas du contrôle, a été montrée d’une fa¬ 
çon explicite par Fleming [16]. Plus tard, ces résultats ont été généralisés 
dans le cas où a est contrôlé par El Karoui [11], 

Remarque : L’injection de l’espace U dans K, est effectuée par 
l’application : 

ip : U —» 3Î 

définie par : 

ip(u)(dt,da ) = dtô Ut (da ) 
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où 5 Ut représente la mesure de Dirac au point u t . 


Exemple illustratif 


Considérons l’exemple du contrôle ordinaire considéré auparavant. 
On identifie u n (t) par la mesure de Dirac S Un(t} sur A . 

Rappelons que u : [0,1] —> A. On considère la mesure q n (dt,du) = 
5 Un (t)(du)dt sur l’espace ([0,1] x A\). 

Lemme 2.8 (Voir Mazliak [30])La suite q n converge faiblement vers 
la mesure 

q(dt ,.) = i[£_i + 5i\{du)dt. 


On peut définir alors un nouveau problème de contrôle associé à 
appelé problème de contrôle relaxé. On considère la dynamique 


<?, 


x 


a _ 


Xq + 



u q (ds, du), 


où le coût associé est donné par J[q) = f^(x^) 2 dt. Notons que les mesures 
sont de la forme q(dt, du) — ô Ut (du).dt. Si de plus q(dt, du) = |[5i (du) + 
û_i (du)]dt, alors on a J(g) = 0 et donc le nouveau problème de contrôle 
admet q comme une solution optimale. 


2.7.1 Position du problème 

On se place sur un espace probabilisé filtré : (O, F, (F t )teR+, P)- 
Soit [0, T] C M + U {+oo} l’ensemble d’indices et U l’ensemble des 
contrôles admissibles à valeurs dans l’espace supposé compact A. On 
considère dans toute la suite le problème de contrôle défini par l’EDS : 

dX t = b{t, X t , u t )dt + cr(i, X t )dB u : t e]0, T] 

X 0 = x, x G M m 

et la fonction coût J définie par J[u) = E\g{X T )] et on garde les mêmes 
hypothèses sur les coefficients qui définissent l’équation d’état et la fonc¬ 
tion coût associée. 

On considère V([0,T) x A) l’ensemble des mesures de probabilité sur 
[0, T] x A tel que la projection sur [0, T] coïncide avec dt, la mesure de 
Lebesgue sur M . l’ensemble V([0,T] x A) muni de la topologie de la 
convergence faible des mesures est un espace métrisable compact. 


Définition 2.9 Toute application mesurable 

fi(.,dt,da) : D—» V{[0, T] x A), 
lu i— » ii[dt, da) = dt, da) 
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est appelée contrôle relaxé. 

L’ensemble des contrôles relaxés est noté par U. 

Proposition 2.10 Etant donné un contrôle relaxé p à valeurs dans 
A. Alors pour tout t G [0, T], il existe une mesure de probabilité p t sur 
A, telle qu’on ait la décomposition suivante : 

da) = dtp(t, da) = dt/j, t (da). (7) 

Preuve : Provient directement du théorème de Fubini.H 

Remarque : La valeur d’un contrôle optimal relaxé est prise d’une 
façon aléatoire sur l’espace A avec la mesure de probabilité /j, t (da), par 
contre dans le cas d’un contrôle optimal ordinaire, on associe la valeur 
u t G A à un instant t G [0, T\. 

Définition 2.11 Soit (p n ) ne ^ une suite de contrôles relaxés. On 
dit qu’elle converge vers fi G ï? (faiblement), si pour tout fonction f 
continue à support compact sur [0, T] x A, on a : 

lim / f(t,a)n n (dt,da) =: / f(t,a)p(dt,da). 

n S -+00 0,T]xA J[0,T\xA 

On constate que la mesure de Lebesgue est la marginale de tous les élé¬ 
ments de ï? sur [0, T], on obtient une convergence stable en affaiblissant 
les hypothèses sur / comme suit : 

Proposition 2.12 Soit ( p n ) une suite convergente vers fi dans 
Alors pour toute fonction mesurable f : [0, T] x A —> R telle que 
Vt G [0,T], l’application 


f{t, •) : A > R 

est continue, on a une convergence stable au sens que : 

lim / f(t,a)p n (dt,da)= / f{t,a)p(dt,da). (8) 
n »+oo 7[ 0r ] xyl J[0,T]xA 

L’ensemble des contrôles relaxés admet une propriété de compacité très 
utile : 

Proposition 2.13(Voir El Karoui [11]) Si on suppose que l’espace 
A est compact, alors ïî l’est aussi. 
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2.7.2 Étude d’un problème de contrôle relaxé 

Dans le problème de contrôle relaxé, on utilise à l’instant t une 
mesure de probabilité sur A (l’ensemble des valeurs de contrôle ) notée 
Ht au beu de u t élément de A. Le problème est alors de modéliser par 
l’équation (9) dans le cas où le terme de diffusion est contrôlé : 

(dX t = ( f A b(t, X t , a)dnt(a))dt + ( f A a(t, X t , a)d/j t (a))dB t , t e]0,T] 
| A 0 = x , x G R m 

(9) 

et par l’équation (10), dans le cas où le terme de diffusion n’est pas 
contrôlé : 

\dX t = (J A b(t,X t ,a)dnt(a))dt + a(t,X t )dB t , t e]0,T] 

| A 0 = x , x G R m 

Etude de problème de contrôle relaxé : a non contrôlé 

On considère les problèmes de contrôle relaxé (10), où le terme de 
diffusion a n’est pas contrôlé. Soient 

b: [0, T] x M m x A —* R m 


et 

a : [0, T] x M m —» At mxd (M) 

deux fonctions boréliennes. Supposons que b est bornée, continue en a 
et uniformément continue en (t, x ) , et que les hypothèses (3) et (4) sont 
vérifiées. 

Sous ces hypothèses, l’équation (10) admet une solution unique forte 
X — (X t ) tGK +, continue, (F t )teR+ adaptée qui X vérifie l’inégalité 

E [| sup \X t \ p \\ < N, Vp> 1, 
te[o,T] 


où N est une constante qui dépend de C, P, T. En outre la solution X 
est définie par : 


Ah 



b(s, X s ,a)/j, s (daŸjds + / a (s, X s )dB s . 


On définit la fonction coût par 


J(H) = E[g(X T )\, 


où Xt est la valeur terminale de la solution de l’équation (10) et g : 
—> R est une fonction de classe C 1 (où C l représente l’ensemble des 
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fonctions continûment différentiables) vérifiant les conditions suivantes : 

i- Il existe une constante positive c, telle que : 

Vx G R m : |g(x)| < c(l + |x|) (11) 

ii- Il existe une constante positive k , telle que : 

\9x\ < k, (12) 


où g x désigne la gradient de g. 

On cherche à trouver un contrôle g* G 3Î minimisant le coût J, autrement 
dit : 

J (fi*) < J (fi), V g G K. 

On note que J(g*) = inf Me sR J(fi). 

Remarques : 


1— Si on pose : 


b(t,x,fi t ) 

l’équation (10) devient : 



b(t,x,a)fi t (da ), 


dX t = b(t, X t , g t )dt + a(t, X t )dB t , t g]0, T] 
X 0 — x, x G M m 


(13) 


et dans ce cas la fonction b satisfait les mêmes conditions que b, et de 
plus elle est linéaire en g . 

On déduit que par l’introduction du contrôle relaxé, on a remplacé l’es¬ 
pace d’état des contrôles ordinaires A par un espace plus large V([0 , T] x 
A) qui est compact et convexe et b par b qui est également linéaire en g. 
2— Comme un cas particulier en mettant g t = d Ut (mesure de Dirac au 
point u t ) pour tout t G [0, T] , on obtient : 


b(t,x t ,fit) = / b(t,x t ,a)a Ut (da) = b(t,X t ,u t ) 

J A 

et on aura un problème de contrôle ordinaire. 

2.7.3 Approximation des trajectoires 

Le problème de contrôle relaxé est réellement une généralisation du 
problème de contrôle ordinaire par le fait que l’infimum de la fonction 
coût sur 5?, l’espace des contrôles relaxés est le même que son inhmum 
sur U, l’espace des contrôles admissibles. 
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Lemme 2.14 (Voir Fleming[16]) Soit fi un contrôle relaxé. Alors il 
existe une suite ( u n ) ne ^ de contrôle ordinaires admissibles de U, telle 
que la mesure dt/j%{da) converge vers la mesure dt/i t (da), où l’on posé : 


Ht(da) = S u n(da). 

En d’autres termes, pour toute fonction 

f : [0, T] x R m x V(A) —» R 

(t,X,fJLt) 1 —» f{t,x,nt), 

fonction continue sur [0, T] x V{A) et linéaire en fi t , on a : 

lira f f(s, X s ,n n s )ds = [ f(s,X 3 ,n s )ds, (14) 

n —*- 00 J 0 J 0 

uniformément en t G [0,T], 


Remarque : Si on note par X n la solution de l’équation (13) asso¬ 
ciée à /j, n pour tout n G N, alors X n doit satisfaire l’équation : 

f dXf = bit, Xf, fj$)dt + a(t, Xf)dB t ; t G ]0, T] 

\XS — x, x E M m . 

Lemme 2.15 Soit n un contrôle relaxé et X la solution de l’équa¬ 
tion (13) associée à /a. Alors il existe une suite (u")„ e ^ de contrôles 
admissibles telle que : 


lim E 


sup \X t -Xf\ 2 

te[o,T] 


0 , 


et 


lim J(p. n ) — J {f. t), 

n—>oo 

où in n ) est une suite de contrôles relaxés, telle que : 

li n (dt,da ) = dtô u n(da), 

et X n est la solution de (13) associée à pA pour tout n G N. 


(16) 

(17) 


Preuve : La démonstration se fait en 2 étapes. 
1— On montre d’abord que 


lim E 


sup \X t 
te[o,r] 



0 . 
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Soient X solution de (13) et X n solution de l’équation (15) associée 
respectivement à /i et à fi. On obtient par différence 


X t -X?= / {b(s,X a ,fjL a )ds-b{s,X^fi))ds+ (a(s,X a )-a{s,XÏ))dB a . 


En ajoutant et en retranchant le terme f*b(s,X a ,/i™)ds, on aura : 


Ab - A7 = J [b(s, X s ,fi 3 ) - b(s, X s , fi )J ds 

+f* (î Ks,x s ,fi)-b(s,x:,fi)y s 

+ (cr(s,A s ) — <j(s,Xg)jdB a . 

et en posant 

^n(^) J* ^^(*^1 x a , Us) 5(s, X s , fjj s )^ ds, 

on obtient : 

Ab - X? = a n (t) + J* ( b(s , Ab ,fi) - b(s, Xfi )) ds 

+ (a(s,X s )-a(s,X:)yB a . 

En appliquant l’inégalité : 

(a + b + c) 2 < 3(a 2 + b 2 + c 2 ) pout tout a,b,c> 0, 


on aura 



\X t -X^<3\a n (t)\ 2 + 3 / ib(s,X a ,fi)-b(s,X?,fi))ds 


1 2 


+3 


a{s,X a )-a{s,X?))dB. 


L J o 


Il résulte alors de l’inégalité de Cauchy-Schwartz que 


\X t -X^<3\a n (t)\ 2 + 3T / \b(s,X s ,fi)-b(s,X:,fi)\ 2 ds 


+3 


J (a(s,X a )-a(s,X?)yB a 


1 2 
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et en passant au sup sur [0, T] et en appliquant l’espérance, on conclut,d’après 
l’inégalité de Burkholder-Davis-Gundy, qu’il existe une constante posi¬ 
tive C, telle que : 


E 

sup \X t -X?\ 2 ' 

<3 E 

sup |a; n (t)| 2 


L te[o,r] 


L te[o,r] J 


+3 TE 
+3 CE 


\b(s,X s ,^ s )-b(s,X:,^ds 


<7 


(s,X s ) - a(s,X?)\ ds 


Puisque b{t , x, u ) et a(t, x ) sont lipschitziennes en x, alors : 


E 

sup \X t -X?\ 2 

< 3 E 

sup | a n (t)\ 2 

+3 TKE 

[ \X S -X?\ 2 ds 

+3 CEE 

' f \Xs-x: 


L te[o,T] 


L ie[0,T] 


L J o 


^ J 0 


On pose D = QK sup(T, C ), d’où 


E 

sup \X t — X™\ 2 

<3 E 

sup | a n {t)\ 2 


L te[o,T] 


L te[o,rj J 


+DE 




En appliquant le Lemme de Gronwall, on obtient : 


E 

sup \X t - A7| 2 

< 3 E 

sup \a n {t)\ 2 

exp(DT). 

(18) 


L te[o,T] 


L te[o,T] 




D’après (4) et puisque b(t, x, fx) est continue sur [0, T] x V(A) et linéaire 
en /i, alors d’après (14), on a : 

lim a n (t) = 0 

n —>-oo 

uniformément en t G [0, T] . Comme b est bornée, alors d’après le théo¬ 
rème de convergence dominée de Lebesgue, on conclut que 

lim E[\a n (t)\ 2 ] = 0, 


et donc 

lim E[ sup |a n (t)| 2 l = 0. 

n >oo ig[o,T] 


Il résulte de (18) que : 


lim E 


sup \X t 

te[o,T] 



0 . 
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2— On va prouver maintenant que : 

lim J(g n ) = J(/i). 

n — >oo 

Soient les fonctions coût associées respectivement à g et à /U n , n G N, 
suivantes : 

J(g) = E[g(X T )\ 
et 

J(ji n ) = E[g(X£)\, 

où Xt et X% sont les valeurs des solutions respectives des équations (13) 
et (15), prises au temps terminal T. 

On a 


I JM - •%") I 2 = \m*T)\ - -E[s(^ï)]| 2 = I E\g(X T ) - 9 (XÎ)]| 2 , 
D’après l’inégalité de Cauchy-Schwartz, on aura : 

\J( P )- J(jn \ 2 <\E\g(x T )-a{x^)f. 


Comme g est dérivable et à dérivés bornés, alors elle vérifie la condition 
Lipschitz avec la constante k, d’où d’après (12) 


\J(g) - J(fjL n )| 2 < | E\g(X T ) - g(XÏ)}\ 2 < kE[\X T - 


< kE 


sup \X t - X 1 

L te[o,T] 


n |2 

t I 



et donc 


0 < | J (fi) - J{g n )\ 2 < kE 


sup \X t 
te[o,T] 



D’après (16) et par passage à la limite lorsque n —> oo, on déduit que 


lim J(g n ) = J (g), 


d’où le résultat désiré.■ 

2.7.4 Existence d’un contrôle optimal relaxé 

On s’intéresse à prouver l’existence d’un contrôle optimal relaxé. 
Dans le cas où le coefficient de diffusion n’est pas contrôlé, on va appli¬ 
quer le théorème de Skorokhod, où l’approche est dû à Fleming. D’autre 
part, on va s’intéresser à l’idée des problèmes de martingale dû à El 
Karoui, pour démontrer l’existence du contrôle optimal dans le cas où le 
coefficient de diffusion est contrôlé. 
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a- Existence d’un contrôle optimal dans le cas où le terme 
de diffusion n’est pas contrôlé 

On constate qu’il existe toujours un contrôle relaxé optimal qui mi¬ 
nimise le coût J(u) sur l’espace fî, on se base sur le théorème de Skoro- 
khod. 

Théorème 2.16 (Théorème de Skorokhod) : Soit (. X n ) n une suite de 
variables aléatoires à valeurs dans telle que Pj^ n la loi de X n converge 

étroitement vers p. Alors il existe un espace probabilisé ^f2, F, et une 

suite de variables aléatoires (X n ) n et une variable aléatoire X définies 
sur cet espace, telle que : 

i — Pjjr = / 1. 

ü - p x n = p xn pour n = 1 , 2 ,... 

iii — X n converge en probabilité vers X par rapport à P. 

Preuve : voir [24] ■ 

Remarque : Puisque X n converge en probabilité vers X, alors on 
peut extraire une sous suite ( X nk ) k telle que : 

lim X nk = X V p.s. 

k—>oo 

Le lemme de Kolmogorov suivant sert à affirmer l’existence d’une ver¬ 
sion continue des processus. 

Lemme 2.17 (Lemme de Kolmogorov) : Soit ( X n ) n une suite de pro¬ 
cessus stochastique continus, définis sur un espace probabilisé (Ll,F,P) 
et à valeurs dans telle que : 

i— il existe deux constantes positives M et 7 telles que : 

W1 < M pour tout ne N 

ii— il existe des constantes a, fi et M k (pour k = 1,2,...), telles que : 

E[\Xfi—X™\ a ] < M k \t— s| 1+/3 pour tout n G N et pour tout t, s G [0, k]. 

Alors, il existe un espace probabilisé F, , une sous-suite de 

processus stochastiques (X„ fc ) et un processus X, qui sont tous conti¬ 
nus et définis sur cet espace, tels que : 

i— Les lois de X nk et X nk coïncident. 

ii— X™ k converge vers X, uniformément sur [0, T] P p.s. 
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En utilisant le théorème de Skorokhod et le lemme de Kolmogorov, 
on énonce un résultat principal concernant l'existence d’un contrôle op¬ 
timal relaxé, qui est donné par le théorème suivant : 

Théorème 2.18 II existe un contrôle relaxé fi* G 3Î tel que : 

J(F*) = inf J(fi) = inf J(u) 

HGvi uÇlU 

Preuve : Montrons d’abord que J(fi*) = inf J(fî). On pose 

Aies? 

1 = in f (19) 

Et soit ( fi n ,X n , B n ) pour n— 1, 2,... une famille telle que 
lim J(fi n ) = ü m E[g(X^)] = l, 

n—>oo n —>-oo 

où X n est la solution de l’équation (15). La famille (fi n , X n , B n ) est 
tendue dans V(A) x C(M + ,M m ) x (7(IR + ,IR d ), d’où d’après le lemme 
de Kolmogorov il existe un espace probabilisé (f2, jF, P) et une suite 
(yü n , X n , B n ) et (p,,X, B) définis sur cet espace tel que : 

i— Pour tout n G N, les lois de (fi n \X n ,B n ) et de (jl n ,X n ,ë n ) coïn¬ 
cident. 

ii— Il existe une sous suite (pA k , X nk , B nk ), telle que : 

lim (^,î nt ,B nfc ) = (/i,X,5) Pp.s. (20) 

72 fc- >-00 

dans l’espace T (A) x C'(M + ,M m ) x C'(E + ,M d ). 

A partir de (i), on déduit que : 

È[g(K')} =s[ 9 (v?)]. 

où E représente l’espérance par rapport à la probabilité P. 

Et par le fait que lim E\g(X^ k )] = l, on aura alors 

n k —>oo L v 

lim Ë[g{Xp)\ =1. 

k —>-oo 

D’autre part, d’après (ii) et la continuité de g on a : 

JÙtPPv?)] =-é[ 9 (à' t )], 

d’où d’après l’unicité de la limite, £ , [g , (Xr)] = l- 
D’après (19), on a : 

J(/2) = inf J (fi) (21) 
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fi est un contrôle optimal et il suffit de prendre fi = fi*. 

Il reste à prouver que X est la solution de l’équation (13). 

En effet : 

X7 k = Z+ f [ b(s, X™ k ,a)fi™ k (da)ds + f [ a(s, X^dB^. 
J 0 JA J 0 J A 

b et a étant des fonctions continues, alors d’après (20) on a : 


lim 

k —>-oo 



0 J A 


b(s, X™ k , a)ii™ k (da)ds — / b(s, X s ,u s )ds P p.s. 


et : 


rt _ __ 

lim / a(s, X™ k ^dB™ k = lim / a(^s,X s ^dB s P p.s. 
k n k >oc 70 


k—>oo 


par conséquent, on aura : 

lim Xfi k — Z + f f b(s, X S1 a)fi s (da)ds + f cr(s,X s )d5 s P p.s. 


k—>oo 


0 JA 


Comme lim X nk — X P p.s., alors : 

k—>oo 


X t = Z+ / b(s,X s ,u s )ds + / cr(s,X s )dB s , 

J 0 J 0 

qui est la forme intégrale de l’équation (13), par rapport au mouvement 
brownien B. 

L’égalité : 

inf J(/i) = inf J(u) 

lie K ueu 

est une une conséquence directe du lemme (égalitéslG, 17). 

b- Existence d’un contrôle optimal dans le cas où le terme 
de diffusion est contrôlé 

Lorsqu’on étudie un problème de contrôle, on a besoin d’estimer le 
coût J(u), qui est une fonction des processus (u t ) et (X t ). C’est pour 
cette raison que la distribution du processus (u t , X t ) est importante , 
d’où la notion de solution faible d’une EDS. 


Solution faible d’une EDS et problème de martingale 

Considérons une EDS contrôlée de la forme (2), comme suit : 

f dX t = b(t, X t , u t ) + a(t, X t , u t )dB t ; t e]0, T] 

| Xo — x, x G M m 


( 22 ) 
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où 

b : [0, T] x R m x A —» R m 
et 

a : [0, T] x W 71 x A —» -M mxd (M) 

sont deux fonctions mesurables continues en x et en a. 

L’objectif du problème est de trouver un contrôle qui minimise une 
fonction coût donnée, comme suit : 


J(u) = E / l(s,X a ,u s )ds + g(XÏ) 


(23) 


où l : [0,T] x x M. —» K et g : R m —» K sont des fonctions 
mesurables. On remarque que ce cas contient le cas déterministe (a = 0). 
Un problème similaire vu précédemment qui n’a pas nécessairement de 
solution, raison pour laquelle on introduit les contrôles relaxés qui sont 
des contrôles à valeurs mesures et l’équation d’état (22) sera de la forme : 


X* — x + 


b(s,a,X s )/j, s (daŸjds + 


a(s,a,X s )/d s (da)^jdB s 

(24) 


Le problème se pose avec le second terme du deuxième membre de 
l’équation (24) car il est une intégrale stochastique. Supposons alors 
que fi n (dt,da ) converge vers n(dt,da), on aura l’estimation suivante : 


(a) : =E[( 
'( " 


! 0 J A 


cr(t,a,X s )n t (da)dB t - j j a(s,a,X t )fj^(da)dB t ^ 


= E 


a 


(t,a,X t )(nt(da) - Ht(da)))dB t 


1 2 


En appliquant l’isométrie d’Itô, on obtient 


(a) = E 


a(t, a, X t )(n t (da) - ^(da)) 


0 L J A 


dt ). 


(25) 


La formule (25) n’est pas satisfaisante. 

Il est intéressant de reformuler le problème sans aucune utilisation de 
l’intégrale stochastique, et ceci va être le rôle principal du problème de 
martingale. 

Solution faible d’une EDS 


Rappelons que la solution faible d’une EDS est définie comme suit : 
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Définition 2.19 Une solution faible de l’EDS 


I dX t = b(t, X t )dt + a(t, X t )dB t , t g]0, T] 
[ X 0 = x G M m 

est une collection d’éléments (Ut, F, ( F t ),P, où : 

i — (f2, F, (F t ), P) est une espace probabilisé filtré. 

ii — (B t ) t > o est un (Ff)—mouvement brownien sur f2. 

iii — (X t ) t >Q est une solution forte de l’équation : 

f dX t = b(t, X t )dt + a(t, X t )dB t , t g]0, T] 

| X 0 = x G MF 1 


(26) 


(27) 


Remarque : L’espace de probabilité et le mouvement brownien ne 
sont pas décrits dans le cas d’une solution faible, la solution est cependant 
tout élément (Ul,F, ( F t ),P , (B t ), (X t )). Dans la formulation précédente, 
on ne sait pas comment trouver l’espace et le mouvement brownien 
(B t )t>o défini sur H. Ainsi, il vaut mieux définir une formulation qui 
nous permet d’expliciter convenablement cette distribution, et ceci sera 
fait par les problèmes de martingales. 

Introduisons le générateur associé aux coefficients b et a, qui est défini 
comme suit : 

df df 1 d 2 f 

Lf(t,x) = ^(M) +b(t,x)-^(t,x) + -a{t,x)a\t,x)-^(t,x), (28) 

pour toute fonction / de classe C^’ 2 ([0, T] x R m ). 

Il est facile de vérifier par la formule d’Itô que si (Q, F , (F t ),P, (B t ), (X t )) 
est une solution faible de (26), alors pour toute fonction / G C^’ 2 ([ 0, T] x 
M m ), alors le processus 

f{t,Xf)~ f Lf(s,X 9 )ds, (29) 

J o 

est une (F t )— martingale. 

Définition 2.20 L ’espace canonique C est l’ensemble des fonctions 
continues de [0 ,T] vers K m . 

Le processus (Xi) teR + défini sur C par : 

Mf) = fit), 


est appelé le processus canonique. 

La filtration canonique est définie par C t = a(X s ,s < t). 

La distribution de ( X t ) t eR+ est une mesure de probabilité sur C. 


45 



Si (f2, F , (F t ), P, ( B t ), (X t )) est une solution faible de (26), alors la distri¬ 
bution Px de (Xt)tëM+ est telle que pour toute fonction / G C£’ 2 ([0, T] x 
MF), le processus défini sur C par f(t,X t ) — f Q Lf(s, X s )ds est une 
(C t ) —martingale sous la probabilité Px- On a juste écrit la propriété de 
martingale autrement (29). 


Définition 2.21 Une mesure de probabilité Q sur C est dite solution 
du problème de martingale associé au générateur L si pour toute fonction 
f G Cl' 2 {{{), T] x R™), le processus défini par 

f(t,X t )~ f Lf {s, X s )ds, 

J o 

est une [Ct) —martingale sous la mesure Q. 

On a vu qu’à chaque solution d’un problème de martingale est associé à 
une solution faible de (26). Le théorème suivant affirme que l’inverse est 
aussi aussi vrai et qu’il existe une bijection entre les solutions faibles et 
les solutions d’un problème de martingale. 


Théorème 2.22 (Voir Ikeda Watanabe[24]) : On suppose que Q 
est une solution du problème de martingale. Alors il existe une solution 
faible (fî,F,(F t ),P,(B t ),(X t )) de (26) telle que (X t ) admet Q comme 
distribution. 

Remarque : La formulation de problème de martingale est intéres¬ 
sante car elle nous permet d’obtenir aisément des résultats limites. 

Proposition 2.23 Supposons que b et a sont des fonctions bornées et 
continues en x. Alors l’ensemble des solutions du problème de martingale 
associé, est un ensemble fermé. 

Preuve : Soit (Q n ) n une suite de solutions du problème de martin¬ 
gale convergente vers Q montre que sous Q, pour tout / G Cl’ 2 ([0,T] x 
M m ) le processus défini par 


f(t,X t )~ f Lf (s, X s )ds 
J o 

est une (Ct)—martingale. 

Soit r < t et soit h une fonction continue bornée et CV—mesurable 
sur C. Alors, on a pour tout entier n : 

Qn 

(30) 

Comme la fonction (f(r,X r ) — Lf(s, X s )ds\h, est bornée sur C et 
continue en r sur [0, T], alors on a d’après la définition de la convergence 


(/(*,**)- f Lf(s, X s )ds^ h =Q n ^f(r,X r )— j Lf(s,X s )ds^jh 
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faible de (Q n ) vers, Q lorsque n 


oo et l’égalité (30) tend à : 


Q 




Lf(s, X s )ds^ h 


Q 


ce qui achève la démonstration.■ 



Lf(s,X s )ds s jh 


Remarque : On peut utiliser la compactification dans les problèmes 
de contrôle stochastique relaxé pour assurer l’existence d’un contrôle op¬ 
timal dans le cas où le terme de diffusion est contrôlé. 
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Chapitre 03 

G-contrôle stochastique optimal 

strict 
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3 G-contrôle stochastique optimal strict 

3.1 Motivation 

On considère une équation différentielle stochastique contrôlée, gou¬ 
vernée par un G—mouvement brownien. Le problème de contrôle optimal 
strict a été étudié, en 2014, par Biagini et al [6] ,où ils ont considéré le 
problème de contrôle stochastique sous la condition de volatilité incer¬ 
taine selon l’approche de l’espérance sous linéaire et celle de l’analyse 
quasi-sûre. En appliquant le principe du maximum, leur objectif était 
de trouver les conditions nécessaires et suffisantes d’optimalité vérifiées 
par le contrôle optimal strict en supposant qu’il existe. Sous des condi¬ 
tions de convexité des coefficients, ils ont maximisé l’hamiltonien, ce qui 
revient à optimiser le critère de performance J. 

3.2 Le principe du maximum classique 

Le principe du maximum représente un outil fondamental dans la 
théorie du contrôle, il sert à trouver la commande optimale permettant 
d’amener un système dynamique d’un état à un autre, en présence de 
contraintes portant sur l’état ou les commandes d’entrée. Le principe 
examine la minimisation d’un hamiltonien sur U, l’espace des contrôles 
admissibles. Historiquement le travail initial dans ce domaine est dû 
à Pontryagin, il était centré sur la maximisation d’une fonctionnelle de 
bénéfice, dans le cas déterministe, plutôt que la minimisation d’une fonc¬ 
tionnelle de coût. 

L’objectif du principe de Pontryagin dans le cas classique est d’éta¬ 
blir les conditions nécessaires d’optimalité alors que les équations de 
Hamilton-Jacobi-Bellmann, notées HJB permettent d’énoncer les condi¬ 
tions suffisantes d’optimalité. 

On considère l’EDS contrôlée : 

f dX t = b(t, X tl u t )dt + a(t, X t )dB t , t e]0, T] 

= x E R m 

On cherche à déterminer les conditions nécessaires d’optimalité vé¬ 
rifiées par une solution forte d’une G-équation contrôlée. 

En tenant compte de l’existence d’un contrôle minimisant le critère 
de performance J : 


r T 


J{u)=E / l{t,X t ,u t )dt +g{X T ) 


o 
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il reste à déterminer certaines conditions satisfaites par le contrôle opti¬ 
mal dites nécessaires. Partant de la condition d’optimalité d’une fonction 
réelle différentiable / : K —» R, si x* G M est un minimum de /, alors 
fd(x*) ^ 0, où fd représente la dérivée de /. On déduit que x* est aussi 
solution du problème de maximisation de H(y) donnée par : 

H(y) = -fd{x*)f{y),Vy > x. 

Généralisons cette condition au problème de contrôle, on établit la 
différentielle de la fonction coût J. On suppose alors que J atteint son 
minimum pour un certain contrôle û G U, alors û doit satisfaire les 
conditions nécessaires d’optimalité nommées principe du maximum sto¬ 
chastique. 

Notons que l’utilisation de ce principe est difficile dans le cas stochas¬ 
tique car le processus adjoint qu’on va définir dans la suite est représenté 
par une moyenne conditionnelle. La méthode utilisée dans cette partie 
(dûe à Kushner) est basée sur le fait que la dérive et le coefficient de dif¬ 
fusion sont différentiables par rapport à x pour tout (t,a) G [0, T] x A. 
Supposons que û est un contrôle optimal, alors il existe un processus 
(pt)teK+> tel que pour tout t G [0, T] et pour tout u G U on a la condi¬ 
tion nécessaire suivante : 

H(p h Xf,û)^H(p t ,Xf,u), 

X u étant la solution de l’EDS contrôlée associée au contrôle optimal û , 
H est le Hamiltonien du système et {pt)te k+ est le processus adjoint. 

3.3 Le principe du maximum pour une G— EDS 


On se place dans les conditions du premier chapitre, où l’on considère 
un G—mouvement brownien réel (B t ) t>Q défini sur un espace d’espérance 

non linéaire ^f2, Ti, E^J. 

3.3.1 Position du problème 


On considère le problème de contrôle défini par la fonctionnelle 

T 

J {u) := Ê[ f f(t,X u (t),u(t))dt + g(X u (T ))] = sup J p («) (31) 

J PeP 


où l’on a posé 


T 

J p (u) = E p [ J f(t,X u (t),u(t))dt + g(X u (T))]. 

o 
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L’objectif est de trouver un contrôle optimal réalisant 

J(û) = sup J(u ) 

uEU 

On considère les contrôles u à valeurs dans un ensemble ici. Soit 
X(t) = X u (t ) un processus contrôlé défini par : 


dX(t) = b(t,X(t), u(t))dt + g(t, X(t ), u(t))d (B) t + cr(t, X(t), u(t))dB(t ) 
I(0) = iCl pour 0 < t < T 

(32) 

On suppose que les coefficients a, b, p sont lipschitziens et continus 
uniformément en ( t,u ). De plus si ces coefficients ne sont pas détermi¬ 
nistes, ils doivent appartenir à Mq(0,T) pour tout (x,u) E R x A. 

Soient / : [0, T] xRxil —* M et g : M — > K, deux fonctions mesurables 
telles que / est continue par rapport à x et g est bornée, différentiable 
vérifiant la condition suivante : il existe une constante C > 0 et e > 0 


vérifiant l’inégalité 

\g\x)\ < (7(1 + |x|) 1 + e / 2 

On note par U l’ensemble des contrôles admissibles et on suppose 
que pour tout u E U, la condition d’intégrabilité suivante est satisfaite 

< oo. 


E 


f(t,X(t),u(t))dt 


La fonctionnelle de performance correspondante à une mesure de 
probabilité P EV est définie par 


J p (u) = E p 


' T 

J f(t,X(t),u(t))dt + g(X(T)) 

_o 


(33) 


Pour étudier le problème de contrôle stochastique optimal, on s’inté¬ 
resse à trouver un contrôle û qui vérifie : 

sup J v (u) = J v {û) VP G V (34) 

uEU 

On définit l’hamiltonien par la formule suivante : 


H{t,x,u,p,q ) = f{t,x,u) + 


b(t , x , u) + fi(t, x , u 


d(B) t 

dt 


, d (B). 
p+a{t, x, u) q 
dt 
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et la G—équation différentielle stochastique rétrograde (EDSR) associée 
aux processus adjoints p(t), q(t ), K(t) par 


dp(t) = — —(t)dt + q(t)dB(t ) + dK(t ); 0 < t < T 
ox 

p(T) = g'(X(T)) 


(35) 


Rappelons que la solution de cette G—EDSR existe vu les hypothèses 
sur / et g et par la définition des contrôles admissibles. 

Théorème 3.1 (Biagini et al [8]) : Soit û G U avec la solution cor¬ 
respondante X(t), p(t), q(t), K(t ) relative à l’équation (32). On suppose 
que K = 0 dans l’équation( 35), les fonctions (x,u) —> H(t,x,u,p(t),q(t )) 
et x —» g(x) sont concaves pour tout t, presque sûrement et 


E 


q{t)) \ u=ü( t) | 


Tu- 


{t-&y 


= 0 q.s. 


pour tout t. Alors u = û est un contrôle optimal pour le problème (34). 

3.3.2 Principe du maximum pour le cas d’information com¬ 


plète 

Les conditions de concavité ne sont pas satisfaites dans la majorité 
des applications, le principe du maximum ne nécessite pas cette condi¬ 
tion. De plus la notion du fait que K soit une G—martingale disparait 
de Téquation adjointe. Biagini et al [6] ont prouvé un résultat qui ne 
dépend pas de la concavité de l’hamiltonien. On donne dans ce chapitre 
leurs principaux résultats sans démonstrations. On suppose dans ce qui 
suit, que les hypothèses suivantes sont satisfaites : 

Al. Pour tout u, P G U , (3 bornée, il existe ô > 0 tels que u+a/3 G U , 
pour tout a G (— ô, ô). 

A2. Pour tout t, h avec 0<i<t + /i<Tet pour toute variable 
aléatoire a G L^(D t ), le contrôle /3(s) := al[ t)t+ /j](s) appartient à U. 

A3. Pour tout u, fd G U tels que soit /3 borné, le processus dérivé 
Y(t) défini par 

/ Y(t) := ^-X u+a Ht) 

\ doc •> 

\ Y( 0) = 0 
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existe, où Y est solution de la G —EDS suivante : 

dY(t) = {Ë (i)r(i) + 1 V t)m } dt + { + 


Lemme 3.2 (Voir Biagini et al [6]) : Soit û un contrôle optimal 
associé à la fonctionnelle J v pour une mesure de probabilité P E V et 
soit l’équation adjointe définie par l’EDSR sous P : 

B H 

dp F (t ) = - —(t,X(t),p F (t), <f{t))dt + <f(t)dB{t ); 0 < t < T 
p ¥ (T) = g\X(T)) P—a.s 

dH p 

Alors on a — (t) \ u =û(t) — 0, où Ton posé 
ou 

H p (t):=H(t,X,u,f(t),m) 

Théorème 3.3 (Biagini et al [6]) : Soit û un contrôle optimal associé 
à la fonctionnelle J F pour toute mesure de probabilité P E V. Alors 
l’équation adjointe est donnée par la G—EDSR suivante : 

dfP(t) = —— (t, JV (t), (t), (t))<it + q^itfdBit) + dKitfi 0 < t < T (36) 

ox 

pF(T) = g\X(T )) q.s. 

Théorème 3.4 (Biagini et al [6 \).On suppose que : 

1) 6 = 0, p,(t,x,u) = xum(t) et a(t,x,u ) = xus(t), où les processus 
m et s sont bornées pour t G [0, T], quasi-continus. De plus la capacité 
de l’ensemble {.s( t) = 0} est nulle pour tout t G [0, T], 

2 ) / = 0 . 

Soit û un contrôle optimal pour la fonctionnelle J v pour toute mesure 
de probabilité P G P. Alors on a : 

pj fjG a ^ 

——(t) := —Hfit.X^^u^fit)^^)) = 0 q.s. (37) 
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Chapitre 04 

G-contrôle stochastique optimal 

relaxé 
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4 G—contrôle stochastique optimal relaxé 

Le but de ce chapitre, qui est l’essentiel de cette thèse, est d’étu¬ 
dier le contrôle optimal des systèmes soumis à l’incertitude ou à l’am¬ 
biguïté du modèle en raison d’informations incomplètes ou inexactes, 
ou de concepts vagues. Le climat ou les conditions météorologiques et 
les marchés financiers sont des domaines typiques où l’information est 
soumise à l’incertitude. Par exemple dans les problèmes de choix de por¬ 
tefeuille optimal en finance où les processus de volatilité et de prime de 
risque sont inconnus et difficiles à estimer à partir de données fiables, 
nous devons considérer une famille de différents scénarios. Pour ce faire, 
on doit effectuer une optimisation de portefeuille robuste qui survit à 
tous les scénarios donnés. Les aspects d’ambiguïté telle que la volatilité 
incertaine ont été étudiés par Peng [33,34, 35], qui a introduit un espace 
d’espérance sous linéaire à l’aide du G—mouvement brownien. Par la 
suite, Denis et Martini [9] ont proposé une structure basée sur l’analyse 
quasi sûre à partir de la théorie de potentiel pour construire une struc¬ 
ture similaire en utilisant une famille étroite de mesures de probabilités 
mutuellement singulières. Ces deux approches sont différentes, Denis et 
al [10] donnent une représentation duale pour l’espérance sous linéaire 
avec le G—mouvement brownien comme le suprémum des espérances 
ordinaires sur V. 

Dans ce chapitre, on considère une équation différentielle stochas¬ 
tique contrôlée et on montre l’existence d’un contrôle u G LA à valeurs 
dans un ensemble d’actions U tel que : 

J(u) = inf J{u). 

u€.U 

Ce problème a été étudié par M. Hu et al [21], Biagini et al [6] et A. 
Matoussi et al [29], où les auteurs ont proposé les conditions nécessaires 
et suffisantes d’optimalité, en d’autres termes le principe de Pontryagin 
et la programmation dynamique respectivement. 

En l’absence d’hypothèses de convexité, le problème de contrôle op¬ 
timal peut ne pas avoir une solution parce que U est trop petit pour 
contenir un minimisant. On cherche à trouver un espace 1Z qui contient 
U l’ensemble des contrôles stricts. Pour obtenir ce résultat principal on 
construit, comme dans le cas classique, l’ensemble 7Z des contrôles relaxés 
comme étant un sous ensemble de l’espace des mesures de probabilités 
sur l’espace d’actions U et on montre que 

inf J{u) = inf J{/i) = J {JT), 

uÇjA [iÇzR. 

OÙ 

ji — arg min J {(T). 
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Pour cela, on montrera l’existence d’un contrôle optimal relaxé pour 
une G—EDS. 

4.1 G—espérance et G—mouvement brownien cano¬ 
nique 

Soient H := {w G G(M+,E d ) : cu(0) = 0}, équipé de la topolo¬ 
gie de la convergence uniforme sur les intervalles compacts, B(fl) la 
a —algèbre de Borel associée, fl t := {w . At : w G 11} et B le proces¬ 
sus canonique. Rappelons que B est un G—mouvement brownien. Soit 
F := ¥ b = o la filtration engendrée par B, qui est seulement 

continue à gauche et soit la filtration limite à droite F + := {J- + ,t > 0}, 
définie par Tf := J 7 t + := n s>t J s . 

Pour toute mesure de probabilité P sur (Q, B (fl)), on note par Ef := 
E, + V (Xf ) (resp. Xf := X+V J\f F (Xoo)) la filtration continue à droite 
P—complétée de Jq + (resp. J 7 *,), où J\f v (X+) (resp. A^ p (P 7 00 )) est l’en¬ 
semble des événements P-négligeables de la cr-algèbre Xf (resp. X^). 
Les filtrations respectives seront notées F p et F p . 

Lemme 4.1 (Soner et al [42]) : Soit P une mesure de probabilité 
arbitraire sur (fl, J- T 00 ). Alors on a : 

1) Pour toute variable aléatoire £, ÎFf — mesurable, il existe une Pp.s. 
unique variable aléatoire £, Xt~mesurable telle que £ = £, P p.s.. 

2) Pour tout processus X , F p —progressivement mesurable, il existe 
un processus unique X , F —progressivement mesurable tel que X = X, 
dt x P p.p..De plus, si X est P p.s.continu, alors X peut être choisi P 
p.s. continu. 

Rappelons que la G—espérance, qui a été largement discutée dans le 
premier chapitre, a été définie par Peng [33, 34, 35] à partir de l’équa¬ 
tion de la chaleur non linéaire. Peng [33, 34] construit, à partir d’une 
G—espérance E , une G—espérance E définie sur PL := Lip(M. d ) l’en¬ 
semble des fonctions lipschitziennes sur W 1 , consistante de la manière 
suivante : 

Soit Lip(fl) l’ensemble des variables aléatoires de la forme 
£ = P(-5p, B t2 ,..., B, n ) 

pour une fonction continue bornée lipschitzienne ip sur (M d ) et 0 < 
ti < t 2 < • • • < t n < T. Le processus des cordonnées (B t , t > 0) est un 
G—mouvement brownien et la variable aléatoire B i est G—normalement 
distribuée sous E et pour tout s, t > 0 et G, i 2 , ■ • •, t n G [0, t] on a : 

E[<p(B t i,..., B tn , B t+S — B t )\ = E[ip(B tl , ..., B tn )\, 
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où ip(xi ,..., x n ) = E[cp(xi ,..., x n , \fsBij]. Cette propriété entraine en 
particulier que les accroissements du G—mouvement brownien sont indé¬ 
pendants et que B t+S — B t et B s sont identiquement iV(0, sE)—distribués. 

Pour tout p G [0, +oo), on note par L P G {ÇÏ) la fermeture de Lip{Çi) 
sous la norme 

(bmi)*- 

Lq(Q) ainsi défini est un espace de Banach. Pour tout t > 0, on pose 
Lip{ü t ) := Lip(ü) n L°(ü t ), L p G (ü t ) := L P G {Ü) n L°(ü t ), 
où L°(f2 t ) est l’ensemble des variables aléatoires T t — mesurables. 

4.2 G—intégrales stochastiques 

Pour tout p G [0,+ 00 ), on note M G p (0,T) l’espace des processus 
élémentaires F—progressivement mesurables, à valeurs dans W* de la 
forme 

n—1 

V(t) = Mu.U+ifo), 

i=0 

avec 0 = to < h < • • • < i n -i < t n = T, n > 1 est une subdivision de 
R + et rji G Lip{fl ti ). Soit M G (0,T) la fermeture de M G p (0,T ) sous la 
norme 

m p m p g (o,t) := ^[J mm. 

Pour tout p G M g 2 (0,T), l’intégrale G—stochastique est définie ponc¬ 
tuellement par 


rt N-l 

hi'n) = / Vsd G B s := V Vj(B tMj+1 - B tMj ), 

Jo 3=0 

avec I{rj) := la fonction I : M G 2 (0,T) —>■ L G (Ü T ) est continue et 

peut être prolongée par continuité à M G ( 0, T). 

Le processus de variation quadratique du G—mouvement brownien 
peut être formulé dans L 2 G {Q , t ) comme étant un processus matrice symé¬ 
trique d x d, continu défini par 

{B)f := B t ® B t - 2 f B s ® d G B s (38) 

J 0 

Notons que la diagonale est constituée de processus non décroissants. Ici 
a <S> b désigne la matrice symétrique d x d, pour a, b G R , définie par 
(a <E> b)x = (a, x) b pour x G R d . 
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On définit l’application J : (0, T) —» Ljk(f2r) : 

pT N ~ 1 

J= r, t d(B)f := £ v t m°„ - (3)°). 

3=0 

On notera par Q l’unique prolongement de J à Mq( 0, T) : 

Q(v)-= [ T Vtd(B >?, rj G Mq(0, T). 

J o 

On a l’isométrie suivante (formulée pour le cas d — 1 pour simplifier) : 


Lemme 4.2 (voir [34]) : On suppose d = 1 et on considère p G 
Mq( 0,T). Alors on a 


'/ fT \ T 


r r T i 

( / p(s)d G B s ) 

= Ê 

/ ^( S )d<£>? 

\J 0 / 


J 0 


4.3 La représentation duale de la G-espérance 


Denis et Martini [9] et Denis et al. [10] montrent la représentation 
duale suivante de la G—espérance en termes d’une famille étroite, fai¬ 
blement compacte V de mesures de probabilité mutuellement singulières 
sur (D, £>(D)). Cette dualité exprime la G—espérance par une espérance 
robuste par rapport à V. On se base sur les travaux cités dans [9] et [10] 
pour la construction explicite de V. Soner et al. [41,42] introduisent une 
analyse pour une telle construction et ses conséquences sur la G—analyse 
stochastique et en particulier le problème d’agrégation des processus. 

Proposition 4.3 (voir. [9,10]) Il existe une famille des mesures de 
probabilité faiblement compacte V sur (D,S(D)) telle que pour tout f G 
Lh(V), on a Ê[£] = sup PGP E P [Ç], 

On considère la capacité de Choquet régulière définie dans le premier 
chapitre 


c(A) := supP(4), A G 13(fl). 

PeV 


On note par Afp := Dpe-p-^^oo) les P-ensembles polaires. 
Définition 4.4 Une mesure de probabilité Q sur (fl,B(fl)) est dite 
absolument continue par rapport à V si Q(4) = 0 pour tout A G Afp. 


On considère la célèbre filtration universelle F p définie à l’aide des 
mesures de probabilité mutuellement singulières P G P (voir. Soner et 
al [42]) : 

¥ v := {fFj’} t > o avec fFf Q (jf V Afp) pour t > 0 (39) 

Pev 
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La formulation duale de la G—espérance propose l’aspect d’agrégation 
suivant : 


Lemme 4.5 [Proposition, Soner et al [41]/ 

Soit rj G Mq( 0,T). Alors, p est Itô-intégrable pour tout P G P. De 
plus, pour tout t G [0, T], 

f p(s)dcB s = f p{s)dB s P p.s. pour tout P G V, 

J o J o 

où l’intégrale du second membre de cette égalité est définie au sens 
d’Itô usuel. Par conséquence, le processus de variation quadratique (B) G 
défini par (38) coïncide avec le processus de variation quadratique usuel 

q.s.. 

Dans la suite, on omettra le symbole G de la G—intégrale stochas¬ 
tique et de la G—variation quadratique. 


Parmi les résultats de stabilité étudiés, voici un lemme qui joue un 
rôle important dans notre analyse. 

Lemme 4.6 Si {P n }^i C V converge faiblement vers P G TL Alors 
pour tout £ G Lq(ÇI t ), i? Pn [£] converge vers G p [£]. 


A partir des propriétés du processus de variation quadratique (B) et 
de la formulation duale de la G—espérance, on a les estimations de type 
Burkholder-Davis-Gundy, formulées en une dimension pour la simplifi¬ 
cation de la manière suivante : 

Lemme 4.7 (voir [17]) On suppose que d = 1. Alors pour tout p > 
2 et pour tout p G Mq( 0,T), il existe une constante C p qui dépend 
uniquement de p et T tel que 


- 

PU 

p- 


'/ G \ p/r 

sup 

/ rj r dB r 


< C P Ê 

( / M 2 dr) 


Js 


_\Js J 


<C p \t - s|a 1 / Ê[\p r \ p \)dr. 


Si a est une constante positive telle que d( '^ t < a q.s., alors on a 
pour tout p > 1 et pour tout p G Mq( 0, T), 


sup 

S<1t<£ 


p r d(B) 7 


< a p \t — s\ p 1 / Êfl^l^dr. 


4.4 L’espace des G—contrôles relaxés 

Soit ( U , d) un espace métrique séparable et V(U) l’espace de mesures 
de probabilité sur U muni de sa cr-algèbre de Borel B {U). La classe 
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M([0,T] x U) des G—contrôles relaxés considérée dans ce travail est 
un sous ensemble de M([0, T] x U), l’ensemble des mesures de Radon 
v(dt,da) sur [0, T] x U ayant pour projection sur [0, T] la mesure de 
Lebesgue dt et pour projection sur U une mesure de probabilité pt(da) G 
V(U) (i.e. u(da,dt ) := p t {da)dt). On muni cet espace de la topologie de 
la convergence stable des mesures. La topologie de la convergence stable 
des mesures est la topologie la plus line qui rend l’application 


q 1-4 


<p(t, a)q(dt, da ) 


>u 


continue, pour toute fonction mesurable bornée ip(t , a) continue en a, 
pour tout t fixé. 

Muni de cette topologie, M := M([0, T] x U) est un espace métri- 
sable séparable. De plus, M est compact si U l’est. On notera que la 
topologie de la convergence stable des mesures entraine la topologie de 
la convergence faible (Voir N. El karoui[12,13] pour plus de details). 

On introduit la classe des G—contrôles stochastiques relaxés sur l’es¬ 
pace (fl T , Lip(JÎ T ), E). 


Définition 4.8 Un G—contrôle stochastique relaxé sur (f2y, LipiUlx), E) 
est une mesure aléatoire q{uj,dt,da) = p t (oj,da)dt telle que pour tout 
A G B{U), le processus (pt(A))t£[o,T] est ¥ p —progressivement mesurable. 

On notera que dans ce cas, l’application [0, 1 ]xOg [0,1] définie par 
(s, eu) i-G p s (uj,A ) est £>([0, t]) ® .Ff—mesurable pour tout t G [0 ,T). En 
particulier, le processus (p t (A)) t€ [ 0 ,t] est adapté à la filtration universelle 
¥ v donnée par (39). 


Dans toute la suite, on désigne également par Tl la classe des G—contrôles 
stochastiques relaxés. Comme dans le cas classique, l’ensemble U([ 0, T]) 
des G—contrôles stricts constitué des processus u , ¥ p —adaptés à va¬ 
leurs dans l’ensemble U , s’injecte naturellement dans H via l’application 
$ : U([0, T])—^ 7 1 définie comme suit : 

${u)(dt,da) = ô u (t){da)dt. G 7 1 (40) 

Le lemme suivant est une extension du célèbre Chattering Lemme. En 
fait, ce lemme annonce que chaque G—contrôle relaxé peut être approché 
par une suite de G—contrôles stricts. 

Lemme 4.9[G—Chattering Lemme] Soient ( U,d ) un espace mé¬ 
trique compact séparable et (ft)t >o u n processus ¥ p —progressivement 
mesurable à valeurs dans V(U). Alors il existe une suite (u n (t )) n >o de 
processus ¥ p —progressivement mesurables à valeurs dans U telle que la 
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suite des mesures aléatoires ô Un ( t )(da)dt converge au sens de la conver¬ 
gence stable (donc faiblement) vers fi t (da)dt quasi-sûrement. 

Preuve : On se donne (nt)t> o un G—contrôle relaxé F 73 —progressivement 
mesurable. La construction trajectorielle détaillée de la suite approximée 
(d Un (t){da)dt ) n >o de n t (dà)dt dans N. El karoui[12] s’étend de sorte que 
les G—contrôles stricts soient F p —progressivement mesurables. Le dé¬ 
veloppement de la démonstration est le suivant : 

On considère l’espace (Cl, H, ( Tf V N-p) t>0 ,E ), pour chaque proba¬ 
bilité P G P et on démontre le lemme pour P quelconque. 

Dans le cas où Tf’ = r> e -p( T t P V N-p) i>0 , l’ensemble où le lemme 
n’est pas réalisé est négligeable par rapport à toute probabilité P donc 
il est polaire. Il résulte alors que le lemme est vérifié quasi sûrement. 

Soit Vt l’ensemble des mesures sur [0, T] x U dont la projection sur 
[0, T] coïncide avec la mesure de Lebesgue. Alors Vr est métrisable com¬ 
pact. La densité des mesures de Dirac au points (u 7 )) (contrôles stricts) 
est utilisée dans la démonstration selon la construction suivante : 

ô u n(da ) —> (nt)(da) dt.ô u n(da) —> dt.fi t (da), 

Comme l’espace U des valeurs du contrôle est supposé compact mé¬ 
trique, alors on peut choisir une partition (B^),k = 1,2, ...k(n) de U 
constitués de boréliens dont le diamètre inférieur à Ala partition est dans 
l’espace U. Soit ajb 7 un point arbitraire de Bj. n> . On divise l’intervalle 
[0, n] en sous intervalles de longueur 2~ n avec des bornes (T-"\ T^f), on 
écrit (Tj,Tj + i). Soit ô > 0 et soit le processus à valeurs dans Vt ■ 

{ | / fi(s, B k )ds si t> 5 

| / yu(s, B k )ds si t <8 
o 

Le processus J 3) k(t, ff) = ^k(tj)Vj<t<,t j+1 est Â' p —adapté et satisfait 
la relation : 


Vj : = 1 si te [tj,t j+ 1 ] et tj > ô. 

k 

On divise l’intervalle Tj — [tj,tj+ 1 [ en k(n) intervalles T jk (fj) de 
longueur (tj+\ — tj) avec t un point quelconque dans T r Pour tout n, il 
existe un indice j n tel que : 

à A tj n < 5 + 2 
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Soit /3(n, ô) un point arbitraire dans A. On pose : 


_ / a k si te k = 1,2, k(n) 

^ ' \ P{n,5) si t ^ ]e,n[ 

et (fi) = l] e ,n[(t)dt8 a (t) + l[o )£ ](t)dt5/3( n ,5), où ô est la mesure de Dirac 
associée au contrôle. Notons que ô est J 71 —adapté. On montre mainte¬ 
nant que ô^ n,s \n) converge vers //. On se base sur la convergence des 
fonctions f,g avec / continue à support dans [0 ,T] et g est continue sur 
U, leurs modules de continuité sont notés par wj et w g . 

Pour une fonction f(t,a) continue à support compact sur [0, T] x U, 
on utilise la notation f(t)g(a) au lieu de /(t, a). Dans le Chattering 
lemme classique on a la convergence S u n[da)dt —> g, t [da)dt 
Si T est majoré par n, on a 


J f(t)g(a)dtfj,(t,da) — j f(t)g(a)dtô {n ’ ê) (g)(da) 

[0,T]xt/ [0 ,T]xU 

= j f(t)g(a)dtfi(t,da)-(Jf(t)g(P(n,ô))dt+ J f(t)g(a(t))dt) 

[0,T]xt/ [0,e] ]e,T\ 

-- j f{t)g(a)dtfi{t,da) - j f(t)g(/3(n,ô))dt- j f(t)g(a(t))dt 

[0,T] xt/ [0,e] ]e,T] 

Le module de cette quantité est majoré par : 


2e 


OO WJ II oo 


/(*)E / d(a)d(t,da) - g(a k )fi(t,B k ))dt 


[0,T] 


B k 


( 




m 


I f(t)n(t,B k )dt- J f(t)dt 
V e >i al,* 

Le deuxième terme de l’expression précédente est borné par 
et le troisième terme est égale à : 


g m j ni 


(Al) 


5^0(a*i)( J f(t){n(t,B k )-/j, k (t)}dt + ^2 J f(t)fj, k (t)dt- j f(t)dt ) 

k }e,T] j Tj T iik 

= ^}2g(ak)( j f(t){fi(t,B k ) - V k (t)}dt + ^2g(a k )Y^f(tj) j {dk(f) ~ lj,k(f)} dt 

k j £ ^ T j k j T . 

-^2d{a k )Y^ I {f(t)-f(t j )}dt + '^2g(a k )'^2 j {/(t) - /(t,-)} fi k (t)dt 

k j T jtk k j Tj 


(A3) 
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= _ EE Qfc )E / {/(*) ~ [ U(t) - f(tj)} Vk(t)dt 

k j T jik k i Tj 

La valeur absolue des deux derniers termes est majorée par 2 \\g\\ oc Twf(-A). 
En fait, le premier terme de (A 3 ) vaut : 


f(t) ( n(t,B k ) - - fi(s, B s )ds | dt 


B k ) - - f(s)n(s,B k )ds ] dt 


t-s 


t—S 


T t 

+ / \dt [ fi(s,B k )(f(s) - f(t))ds (A4) 


£ t—S 


On note par C( x le premier terme de ( A 4 ) et Cf" le deuxième terme 
de (A^). D’après le théorème de Fubini, on a 

\C* | < ÔWf f fi(s,B k )ds et |X)fcfl , ( a fc)- C 2 C | < TS \\g\\oo w f 

0 

On peut écrire C^ de la manière suivante : 

T £ 

J f(s)/jL(s,B k )(l-^-j^-)ds+ J f(s)/i(s,B k )^-^ — -ds 

T-S £-5 

on a aussi J2 k g(a k )C^ < 25 H/H^. 

Ainsi, le premier terme de (A 3 ) est majorée par : 

s\\g\LW\L + Tv>f) 


Le deuxième terme de (A 3 ) vaut. 



7 j,k(t))dt = S ^ 2 / g(a k ) E f(tj) / (g k (t) - g k (tj))dt 


t. 


E^E^h / dt / (jJ>(s,B k )~ n(s - ô,B k ))ds 


k j 

= C 3 


T, t. 
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Par une intégration par parties on obtient : 


C 3 = /t(L+i -t)(ji(t,B k ) - n(t - ô, B k ))dt 

k j J 

1 3 

= I J(L+i - t )(529{ak)){lJ‘(t, Bk) - n(t - ô, B k ))dt 

j i k 

Il résulte que 


|C 3 |< n2-^ n \\g\U\f\\ oo /ô (As) 

En sommant les différents résultats : ( A\) — (4 5 ), on obtient : 

v(fg)-à in, 5 \v)(fg)\ < Iklloo |^I Tw / + 2 ( £ + 5 + ) ll/lloo+ Tw /( 5 )} 

+1 WfWiwg 


En prenant ô = ^,ô n = 5"’", on a la convergence vers 0, pour tout 
P G V, d’où dtS Un (t)(da ) converge vers dt/j, t (da) presque sûrement par 
rapport à P. On conclut que dtô Un (t)(da ) converge vers dt/i t (da ) quasi 
quasi-sûrement d’où le résultat espéré. ■ 

4.5 Le G—contrôle stochastique optimal relaxé 

L’objectif de ce paragraphe est d’établir l’existence d’un infimum de 
la fonctionnelle de performance relaxée suivante : 


J{n) = É Lf(t,x fL (t),a)iJ.t(da)dt + h(x ll (T)) 


(41) 


sur l’ensemble 1Z des G—contrôles relaxés pour la G —EDS relaxée sui¬ 
vante : 


dx^(t) — a(t, x y ‘(t))dB t + / b(t,x >l (t),a)n t (da)dt+ / 'y(t,x IM (t),a)iJ, t (da)d(B) t 

Ju Ju 

(42) 

avec la condition initiale x^( 0 ) = x, où 


b:[0 1 T]xR d xU ^ R d , 
a, 7 : [0, T] x R d x U ->• R dxd 
f : [0, T] x R d x U -> K 
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et 


h : R d -7 M 

sont des fonctions déterministes. 

On notera que lorsque p = ô u , u G ld([0, T]), le processus := x u 
résout la G—ED S suivante : 

f dx u (t) — ait , x u (t))dB t + b(t, x u (t),u(t))dt + 7 (t, x u (t ), u(t))d(B) t , 
|x“( 0 ) = x 

(43) 

De plus, selon l’injection décrite précédemment, on peut écrire J(u) — 
J(S u ). Pour cela, on fait les hypothèses suivantes : 

(Hl) Les fonctions b, 7 et a sont continues, bornées et Lipschitziennes 
continues uniformément en ( t,u). 

(H2) Les fonctions f et h sont continues et bornées. 

Il résulte que sous l’hypothèse (Hl), la G—EDS admet une unique 
solution x ^ G Mq( 0, T), pour tout p G Tl, qui satisfait 

E[ sup \x^(t)\ 2 } < 00 (44) 

0 <t<T 

De plus, 

X 11 := [ [ f(t,x fl (t),a)n t (da)dt + h(x^(T)) e L^(D T ) (45) 

J 0 Ju 

Remarque : Les hypothèses (Hl) et (H 2) sont fortes et peuvent être 
affaiblies. On les suppose pour une présentation simple du résultat prin¬ 
cipal. 

On est maintenant en mesure d’énoncer et de montrer le théorème 
suivant, qui constitue le résultat essentiel de notre travail : 

Théorème 4.10 On a 


inf J(u) = inf J (pi) (46) 

ueU[0,T] /Æ n 

qui se réalisent en p G 7Z. 

Rappelons que 

J (h) = sup J ¥ (p) (47) 

PeP 

et que la fonctionnelle de performance relaxée associée à P G P est 
donnée par 


J r {n) = E v P L/(f,x''((),a)^ ( (da)d( + /i(x' i (r)) 


(48) 


Preuve : Pour montrer ce théorème, on utilise le G —Chattering 
Lemme et les résultats de stabilité pour la G—EDS (42). Selon le G—Chattering 
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Lemme, pour tout G—contrôle relaxé p G P, il existe une suite ( u n ) G 
G([0 ,T]) de G—contrôles stricts telle que ô u n^(da)dt converge faible¬ 
ment vers fi t (da)dt quasi-surement. 

On note par x n la solution de la G—EDS associée avec u n (ou ô u n^ ( da )) 
définie par 


dx n (t ) = a(t , x n (t))dB t + b(t , x n (t ), u n (t))dt + q(t, x n (t),u n (t))d(B) t , 
x n {0)=x (49) 

On a les résultats suivants de stabilité pour la G—EDS (42) et la 
performance J p , pour tout P G "P. 

Lemme 4.11 Pour tout P G P, on a 


1) 


lim E Ÿ 

n—ï oo 


sup I x n (t) 

.0 <t<T 



(50) 


et 


2) De plus , 


lim J p (iP) = J P (yu) 


inf J p (ti) = inf J p (/i), 

weM[o,T] mgk. 


et il existe un G—contrôle relaxé jlp G P te/ que 


J p (/ip) = inf J P (yu). 


(51) 

(52) 


(53) 


Preuve : 

l)Remarquons que vu la propriété d’agrégation du Lemme 4.5, les 
G—EDS (42) et (43) deviennent des EDS gouvernées par une martingale 
continue B , sous tout P G P. La preuve de ce résultat se déduit de [1] 
ou [3]. On en donne seulement une esquisse. 

En utilisant les faits que sous P G P, P est une martingale de pro¬ 
cessus de variation quadratique (B) telle que c t est bornée par 

une matrice d déterministe, symétrique d x d définie positive et que aP 
et x n satisfont les EDS : 

cizP(t) = a(t,x^(t))dB t + / (b(t, aP(t), a)+Ct 7 (t, aP(t), a))m t (da)dt, oP(0) = x, 

Ju 

et 

dx n (t ) = a(t,x n (t))dB t +(b(t 1 x n (t),u n (t)+c t 'y(t 1 x n (t),u n (t)))dt, x n (0) = x. 
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En vertu de l’hypothèse (H 1) sur a, b et 7 , la preuve de(50) repose 
sur les inégalités standards de Gronwall et Burkholder-Davis-Gundy, le 
théorème de convergence dominée et la convergence stable de 8 Un u)(da)dt 
vers fi t (dd)dt. De plus, (51) résulte de la convergence stable 8 Un (^{da)dt 
vers / i t (da)dt , qui est garantie par la continuité et les hypothèses (H 2 ) 
sur / et h. 

La preuve de (52)se conduit comme suit. D’après (40), on a J p (u) = 
J r (S u ), d’où inf Me t/[o,T] d v (u) > infJ p (/i). Soit /i G 7 Z arbitraire. 
D’après le G—Chattering Lemme, il existe une suite des G—contrôles 
stricts ( u n ) C U[0 ,T] telle que ô Un ^(da)dt converge faiblement vers 
Ht{da)dt. Il résulte de (51) que 

J p (/r) = lim J Ÿ (u n ) > iuf J p (u ), 

n-ïoo u&A[0,T] 


et comme p est arbitraire, on en déduit que 


inf J P (p) > iuf J p (u). 

I 1 €R ueU[0,T] 


2)La preuve de l’existence d’un G—contrôle optimal relaxé (53), re¬ 
pose sur les faits suivants : Soit (/G) n >o une suite minimisante de inf ;ieK J p (yu) 
tel que infJ p (yu) = lim J p (yU n ). Les contrôles fi n sont des variables 

n—ïo o 

aléatoires dans l’ensemble compact M, d’où d’après le théorème de Pro- 
horov, la famille de distributions associée est relativement compacte dans 
M. L’hypothèse (H 1) entraine que la suite des processus (/G,x A ‘ n ) n >o 
est étroite dans l’espace M x C([0,T],M d ). Il existe alors une sous suite 
{n nk ,x^ k ) k >Q de (// n ,x^) n >o qui converge faiblement vers (yû,x^), d’où 
(44). En utilisant le théorème d’injection de Skorokhod, les hypothèses de 
continuité (H 2) et le théorème de la convergence dominée de Lebesgue, 
on obtient finalement : 


inf J P (yu) = lim J p (/i" fe ) = </ P (/i), 

[iÇzR, k—> 00 

ce qui achève la démonstration.■ 

L’idée de la proposition suivante est de rendre la limite donnée par 
la formule (50) valide sous l’espérance sous linéaire E. 


Proposition 4.12 Supposons que l’hypothèse (H 1) est satisfaite. 
Soient x^ et x n les solutions associées à /a et u n respectivement. Alors 
on a 


lim È 

n—ïoo 


sup I x n (t) 

0 <t<T 


x»(t) I 2 


0 . 


(54) 


Preuve : Posons := sup 0<i<T \x n (t) — x^{t)Ÿ- Il est clair que 

G L^(D r ) pour tout n > 1. S’il existe un ô > 0 tel que .£[£„] > 
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S, n = 1,2,..., alors on peut trouver une probabilité P„ G P tel que 
E p ^ n ] > ô — -, n — 1,2,.... Le fait que P soit faiblement compact 
nous permet d’affirmer qu’on peut extraire une sous suite {Pn fc }£Li qui 
converge faiblement vers P G P, d’où 

lim E ]p % ] = lim lim E Pn *[£ n ) > lim inf E ?n k [£ n J > ô, 

j—ï oo j—>oo k—ï oo fc—>• oo 

ce qui est en contradiction avec le fait que lim E p [( n .} = 0 d’après le 

j-t oo 3 

lemme 4.11.B 

Remarque la méthode de la démonstration de (50) ne peut pas 
s’étendre à notre cas pour montrer (54), car sous l’espérance sous linéaire 
le théorème de convergence dominée et le lemme de Fatou ne sont pas 
valides, par contre les inégalités de Gronwall et Burkholder-Davis-Gundy 
(voir Lemme 4.7 ) restent valables pour les G—EDS et les G—intégrales 
stochastiques browniennes. 

Corollaire 4.13 Supposons que l’hypothèse (H2) est satisfaite. Soient 
J(u n ) et J (pi) les fonctionnelles de performance correspondantes respec¬ 
tivement à u n et p avec dtS u n^(da) convergente faiblement vers dtpL t (da) 
quasi-sûrement. Alors, il existe une sous suite ( u nk ) de ( u n ) telle que 

lim J(u nk ) = J (/a). 

k—ïoo 

Preuve : D’après Denis et al[10], on obtient l’existence d’une sous 
suite ( x nk (t)) k qui converge vers x^(t) quasi-surement, uniformément en 
t. On a, d’après le Lemme 4.11, pour tout P G P, 

lim J p (u nk ) = J p (pc). (55) 

k—> oo 

En fait d’après l’égalité J (pi) — min J(yu), on note que J(u nk ) = 
E[x uUk } et J (pi) = E\x p ], où les variables aléatoires x et X 1 '' appar¬ 
tiennent à Lg(D T ). S’il existe 5 > 0 tel que E[x uUk ] > E[x' J ] + ô, n*, > 
i, l + 1,..., on peut alors trouver une mesure de probabilité P m G V tel 
que 

E^[ x ^ k ]>Ê(x^}+ô--. 

m 

Du fait que V soit faiblement compact, on peut extraire une sous suite 
{Pmfcl/bLi qui converge vers P G P. On a alors 

E P [x >1 ]~ lim E Pmk [x^ 1 ] = lim lim E Pmk [x Unj ] 

k—ïoo k—ïooj—ïoo 

> lim inf, E Pm i [x Un ' 3 ] > lim inf ( E [x >l ] +5-— J = E [x M ] + ô. 

j->c jo ' j-too y ' rnj J 
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Par suite, E p [x M ] > E[x^]+ô, qui est contradictoire avec la définition 
de l’espérance sous linéaire. Par suite 

lim J{u nk ) < J(fi). 

k—ïoo 

D’autre part, en utilisant (55) on a pour tout P G V 
lim J(u nk ) > lim J p (u nk ) = J p (/r), 

k—ï oo k—> oo 

par suite lim J(u nk ) > J(/jl), ce qui achève la démonstration. ■ 

k—ï oo 

Preuve de Théorème 4.10 On a d’après la définition d’un contrôle 
relaxé et le Corollaire ci dessus : 

inf J(u) < inf 
ueU[0,T] n€iz 

D’autre part, on a pour tout u G U[0, T], ô u G 7 Z. De plus, 

J(u ) = J(S u ) > inf J(yu), 



Il résulte que 

inf J(u) > inf J(n), 

ueU[0,T] 

d’où le résultat (46). 

On revient maintenant à la démonstration de l’existence d’un G—contrôle 
optimal relaxé. Le fait que f et h sont continues et bornées pour tout 
u G Tl 


X u -= [ [ f(t,x v (t) 1 a)v t (da)dt + h(x v (T)) el} G (ü T ). 

J 0 Ju 

Par le lemme 4.6, on obtient alors pour tout v G 1Z, 

lim J Pn (u) = J Q (v), (56) 

n—>oo 

On en déduit que la suite {Pn})^ G V converge faiblement vers Q G V. 
On suppose qu’il existe un £ > 0 tel que, pour tout v G 7Z, 

J(u) > inf J(n) + e. 

fi€TZ 

D’après le lemme 4.11, il existe un G—contrôle relaxé fi G 1Z tel que 
fi p = arg rnin /ieK ,7 lp (q) pour tout P G ? et on obtient 

J(v) > sup inf J P (/i) + £ = sup J p (fip) + £. 

p eV ^ e:R - PeP 
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D’autre part, pour tout n > 1, il existe P n G V tel que 

J Vn (v) > J(u ) + 

n 

La suite {P ?i })^i G T 5 devient faiblement compact, on peut extraire une 
sous suite {P^ ljL! qui converge faiblement vers Q G V. Il résulte de 
(56) que, pour tout v G Tl, 

J Q (v) = lim J ¥n i (is) > sup J P (/ûp) + e. 

j-*<x> pgp 

En particulier, pour un v® G 1Z donné, on obtient 

J Q (u Q ) > J Q (v Q ) +£, 

ce qui contredit le fait que e > ().■ 
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Conclusion et perspectives 

Dans cette thèse, on a prouvé l’existence d’un G—contrôle optimal 
relaxé dans le cadre de l’espérance non linéaire, satisfaisant les conditions 
de régularité nécessaires, en s’appuyant sur les résultats de S. Bahlali, 
B. Mezerdi et B. Djehiche sur le problème de contrôle stochastique clas¬ 
sique réalisés en 2006. Partant d’une équation différentielle stochastique 
gouvernée par un G—mouvement brownien ; ces résultats ont été diffi¬ 
ciles à obtenir à cause de la non-linéarité de la G—espérance et de la 
difficulté de définir un espace adéquat avec la nature du G—contrôle re¬ 
laxé. Pour cette raison, on a considéré la filtration universelle proposée 
par Soner et al en 2012. Ajouté à cela, notre principale difficulté réside 
dans le fait que le théorème de skorokhod n’est pas généralisé dans le 
cadre du G—calcul stochastique, d’où l’utilisation des méthodes de com¬ 
pactification traitées par N. El karoui et al en 1987. On a prouvé cette 
existence en se basant sur le G—Chattering lemme, qui est une généra¬ 
lisation du célèbre Chattering lemme classique, obtenu dans le cadre de 
la G—espérance qui a été la clé de nos résultats. 

On a également prouvé, en utilisant l’approximation des trajectoires, 
que d’une part chaque G—contrôle relaxé est une limite d’une suite de 
contrôles stricts bien définis et d’autre part, que chaque processus de 
diffusion hé au G—contrôle relaxé est une limite au sens fort d’une suite 
de processus de diffusion associés aux G—contrôles stricts. 

Ce travail de recherche a été couronné par la démonstration d’un 
théorème essentiel qui montre que l’infimum du coût sur l’espace des 
contrôles stricts est le même sur l’espace des contrôles relaxés, comme 
dans le cas classique. 

Comme perspectives, nous envisageons de traiter des applications en 
finance dans le cas où la volatilité est incertaine, ce qui est assuré par 
l’existence d’un G—contrôle optimal relaxé et éventuellement considérer 
d’autres problèmes en particulier, les G—équations stochastiques rétro¬ 
grades contrôlées ainsi que le problème de contrôle mean-held. 
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